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{Sj . Abstract 

This presentation reviews an approach to nuclear many-body systems based on the sponta- 
neously broken chiral symmetry of low-energy QCD. In the low-energy limit, for energies and 
momenta small compared to a characteristic symmetry breaking scale of order 1 GeV, QCD is 
realized as an effective field theory of Goldstone bosons (pions) coupled to heavy fermionic sources 
(nucleons). Nuclear forces at long and intermediate distance scales result from a systematic hier- 
archy of one- and two-pion exchange processes in combination with Pauli blocking effects in the 
nuclear medium. Short distance dynamics, not resolved at the wavelengths corresponding to typi- 
cal nuclear Fermi momenta, are introduced as contact interactions between nucleons. Apart from 
a set of low-energy constants associated with these contact terms, the parameters of this theory 
are entirely determined by pion properties and low-energy pion-nucleon scattering observables. 
^^O This framework (in-medium chiral perturbation theory) can provide a realistic description of both 

isospin-symmetric nuclear matter and neutron matter, with emphasis on the isospin-dependence 
determined by the underlying chiral NN interaction. The importance of three-body forces is empha- 
sized, and the role of explicit A(1232)-isobar degrees of freedom is investigated in detail. Nuclear 
chiral thermodynamics is developed and a calculation of the nuclear phase diagram is performed. 
This includes a successful description of the first-order phase transition from a nuclear Fermi liquid 
to an interacting Fermi gas and the coexistence of these phases below a critical temperature T c . 
Density functional methods for finite nuclei based on this approach are also discussed. Effective 
interactions, their density dependence and connections to Landau Fermi liquid theory are outlined. 
Finally, the density and temperature dependence of the chiral (quark) condensate is investigated. 

1 Introduction 

Almost eight decades ago, Yukawa's pioneering article [1] introduced the framework for a systematic 
approach to the nucleon-nucleon interaction, based on the exchange of a boson identified later as the 
pion. In the footsteps of Yukawa's original work the next generation of Japanese theorists already made 
impressive efforts to proceed from the long-range part of the interaction to shorter distances between 
nucleons. A cornerstone of these developments was the visionary conceptual design by Taketani et al. 
[2] of an inward-bound hierarchy of scales governing the nucleon-nucleon potential, sketched in Fig. [TJ 
The long distance region I is determined by one-pion exchange. It continues inward to the region II 
of intermediate distances dominated by two-pion exchange. The basic idea was to construct the NN 

1 



O 

in 



o 



% 




N N 



► r \m r 



two-pion 

exchange N N 

7T A one-pion 

exchange 



Figure 1: Schematic picture displaying the hierarchy of scales governing the nucleon-nucleon potential 
(adapted from Taketani [3]). The distance r is given in units of the pion Compton wavelength /i _1 ~ 1.4 
fm. 



potential in regions I and II by explicit calculation of 7r and 2-k exchange processes, whereas the detailed 
behaviour of the interaction in the short distance region III remains unresolved at the low-energy scales 
characteristic of nuclear physics. This short distance part is then given a suitably parametrized form. 
The parameters are fixed by comparison with scattering data. 

Taketani's picture, although at the time of its appearance not at all universally accepted by the 
international community of theorists, turned out to be immensely useful. Today this strategy is the one 
pursued by modern effective field theory approaches. It is amazing how far this program had already 
been developed in the late fifties of the twentieth century. One example is the pioneering calculation of 
the two-pion exchange potential [U [5] (using dispersion relation techniques) and early knowledge [B] of 
the resonant pion-nucleon amplitude which anticipated the A-isobar models of later decades [7] . 

Today's theory of the strong interaction is Quantum Chromodynamics (QCD). There exist two 
limiting situations in which QCD is accessible with "controlled" approximations. At momentum scales 
exceeding several GeV (corresponding to short distances, r < 0.1 fm), QCD is a theory of weakly 
interacting quarks and gluons. At low momentum scales considerably smaller than 1 GeV (corresponding 
to long distances, r > 1 fm), QCD is governed by color confinement and a non-trivial vacuum: the 
ground state of QCD hosts strong condensates of quark-antiquark pairs and gluons. Confinement is 
accompanied by the dynamical (spontaneous) breaking of Chiral Symmetry, a global symmetry of QCD 
that is exact in the limit of massless quarks. Spontaneous chiral symmetry breaking in turn implies 
the existence of pseudoscalar Nambu-Goldstone bosons. For two quark flavours (Nj = 2) with (almost) 
massless u and d quarks, these Goldstone bosons are identified with the isospin triplet of pions. Low- 
energy QCD is thus realized as a (chiral) effective field theory in which these Goldstone bosons are the 
active, light degrees of freedom. In the low-energy, long- wavelength limit, Goldstone bosons interact 
weakly with one another or with any massive hadrons. In this limit perturbative methods are applicable 
and a systematic expansion in powers of a "small" parameter can be performed (chiral perturbation 
theory jSJE]). 

The chiral effective field theory (ChEFT) just sketched provides the interface between QCD and 
nuclear physics. This theory has become the framework for a successful description of the nucleon- 
nucleon interaction [10J HU [12], and it is the starting point of a systematic approach to the nuclear 
many-body problem and its thermodynamics at densities and temperatures well within the confined 



(hadronic) phase of QCD. This review is written to serve as an introduction and state-of-the-art overview 
of this active area of research. It does not aim for completeness, however, and should be read in 
combination with complementary overviews, such as refs. [TU j ITT j [T2]. 

2 Low-Energy QCD and Chiral Symmetry 

Before turning to the more detailed presentation of the ChEFT approach to nuclear many-body systems, 
it is useful and instructive to briefly recall how the special role of the pion emerges in the context of chiral 
symmetry in QCD, through the Nambu-Goldstone mechanism of spontaneous symmetry breaking. In 
other words: how does Yukawa's pion figure in the framework of QCD? 

2.1 Chiral Symmetry and the Pion 

Historically, the foundations for understanding the pion as a Nambu-Goldstone boson [131 E3] were 
initiated in the 1960 's, culminating in the current algebra approaches [15] (combined with the PCAC 
relation for the pion) of the pre-QCD era. A most inspiring work from this early period is the one 
by Nambu and Jona-Lasinio [2] (NJL). Just as the BCS theory provided an understanding of the 
basic mechanism behind superconductivity, the NJL model helped clarify the dynamics that drives 
spontaneous chiral symmetry breaking and the formation of pions as Goldstone bosons. 

Consider as a starting point the color current of quarks, J^ = qj^q, where t a (a = 1,...,8) are 
the generators of the SU(3) C color gauge group and q denotes the quark fields with 4N c Nf components 
representing their spin, color and flavor degrees freedom (N c = 3). This current couples to the gluon 
fields. Assume that the distance over which color propagates is restricted to a short correlation length 
l c . Then the interaction experienced by low-momentum quarks can be schematically viewed as a local 
coupling between their color currents: 

C mt = -G c Jl(x)JZ(x), (1) 

where G c ~ g 2 I 2 is an effective coupling strength of dimension length 2 which encodes the QCD coupling, 
g, averaged over the relevant distance scales, in combination with the squared correlation length, I 2 . 

Now adopt the local interaction Eg.([T|) and write the following model Lagrangian for the quark fields 
q(x): 

C = q(x)(iYd^ - m q )q{x) + C mt (q, q) . (2) 

In essence, by "integrating out" gluon degrees of freedom and absorbing them in the four-fermion 
interaction £ int , the local SU(3) C gauge symmetry of QCD is now replaced by a global one only. 
Confinement is obviously lost, but all other symmetries of QCD are maintained. The mass matrix m q 
incorporates small "bare" quark masses. In the limit m q —¥ 0, the Lagrangian (J2J) has a chiral symmetry 
of left- and right-handed quarks, SU(Nf) L x SU(Nf) R , that it shares with the original QCD Lagrangian 
for Nf massless quark flavors. 

A Fierz transform of the color current- current interaction (CQ) produces a set of exchange terms 
acting in quark-antiquark channels. For the Nf = 2 case: 
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with the isospin Pauli matrices r — {ti,t 2 ,t 3 ). For brevity we have not shown a series of terms with 
combinations of vector and axial vector currents, both in color singlet and color octet channels. The 
constant G is proportional to the color coupling strength G c . The ratio of these two constants is uniquely 
determined by the number of colors and flavors, N c and Nf. 



The steps just outlined can be viewed as a contemporary way of introducing the time-honored NJL 
model [13]. This model has been further developed and applied [THJ [T7] to a variety of problems in 
hadron physics. The virtue of this schematic model is its simplicity in illustrating the basic mechanism 
of spontaneous chiral symmetry breaking, as follows. In the mean-field (Hartree) approximation the 
equation of motion derived from the Lagrangian (J2J) leads to a gap equation 

M q =m q - G(0|gg|0> , (4) 

which links the dynamical generation of a constituent quark mass M q to the appearance of the chiral 
quark condensate 

(0|w|0> = -tr lim(0|Tg(0)g(x)|0) = -2iN f N e J ( ^ 4 ^f}^ -^ . (5) 

This condensate plays the role of an order parameter of spontaneous chiral symmetry breaking. Starting 
from m q = a non-zero quark mass M q develops dynamically, together with a non-vanishing chiral 
condensate (0|gg|0), once G exceeds a critical value of order G cr ;t — 10 GeV -2 . The procedure requires 
a momentum cutoff A ~ 1M q beyond which the interaction is "turned off" . Note that the strong non- 
perturbative interactions, by polarizing the vacuum and turning it into a condensate of quark- ant iquark 
pairs, transmute an initially pointlike quark with its small bare mass m q into a massive quasi-particle 
with a size of order (2M q )^ 1 . 

2.2 The Pseudoscalar Meson Spectrum 

The NJL model demonstrates lucidly the appearance of chiral Nambu-Goldstone bosons. Solving Bethe- 
Salpeter equations in the color singlet quark-antiquark channels generates the lightest mesons as quark- 
antiquark excitations of the correlated QCD ground state with its condensate structure. Several such 
calculations have been performed in the past with Nf = 3 quark flavors [HIj \TT\ [T5] . Such a model 
has an unwanted U(3)l x U(3)r symmetry to start with, but due to the axial U(1)a anomaly of QCD 
this symmetrjo should be reduced to SU(3)l X SU(3)r x U(1)v- In the NJL model, instanton driven 
interactions are incorporated in the form of a flavor determinant [19] det[^(l ±75)^]. This interaction 
involves all three flavors u, d, s simultaneously in a genuine three-body (contact) term. 

The symmetry breaking pattern resulting from such a calculation is apparent in the pseudoscalar 
meson spectrum of Fig. [2j Starting from massless u, d and s quarks, the pseudoscalar octet emerges as 
a set of massless Goldstone bosons of spontaneously broken SU(3)i x SU(3)r, while the broken U(1)a 
drives the singlet 770 away from the Goldstone boson sector. Finite quark masses shift the pseudoscalar 
(J 71 = 0~) nonet into its empirically observed position, including ^o-^s mixing. 

The very special nature of the pion as a Nambu-Goldstone boson is manifest in the famous Gell- 
Mann-Oakes-Renner relation [20] derived from current algebra and PCAC: 

mlrt = ~(m u + m d )(0\qq\0) + O(ml d \nm u , d ). (6) 

It involves the pion decay constant, f n = 92.4 MeV, defined by the matrix element which connects the 
pion-state with the QCD vacuum via the isovector axial vector current, A^ = 97^75^- q : 

(o\A^o)\Mp)} = ^ abP fl U. (7) 

Just like the chiral condensate (0|gg|0), the pion decay constant /„. is a measure of spontaneous chiral 
symmetry breaking expressed in terms of a characteristic scale A x ~ 4^/^ ~ 1 GeV. The non-zero pion 



lr The axial [/(^-transformations of quarks, q — > exp(ia7s)g, constitute only a symmetry of classical chromodynamics, 
but not of its quantum version. 
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Figure 2: Symmetry breaking pattern in the pseudoscalar meson nonet calculated in the three-flavor 
NJL model 
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138 MeV ^C A x , reflects the explicit symmetry breaking by the small quark masses, with 
q . One should note that the quark masses m u ^ and the condensate (0|gg|0) are both scale 
dependent quantities. Only their product is scale independent, i.e. invariant under the renormalization 
group. At a renormalization scale of about 2 GeV, a typical average quark mass \{m u + m d ) ~ 3.5 
MeV implies (0|gg|0) ~ -(0.36 GeV) 3 . 



2.3 Scales and Symmetry Breaking Patterns 

The quark masses are the only parameters that set primary scales in QCD. Their classification into 
sectors of "light" and "heavy" quarks determines very different physics phenomena. While the heavy 
quarks (i.e. the t, b and - within limits - the c quarks) offer a natural "small parameter" in terms of 
their reciprocal masses, such that non-relativistic approximations (expansions of observables in powers 
of l/mt,b,c) tend to work increasingly well with increasing quark mass, the sector of the light quarks 
(i.e. the u, d quarks and - to some extent - the s quark) is governed by quite different principles and 
rules. Evidently, the quark masses themselves are now "small parameters", to be compared with a 
characteristic "large" scale of dynamical origin. This large scale is visible as a characteristic mass gap 
of about 1 GeV which separates the QCD vacuum from almost all of its excitations, with the exception 
of the pseudoscalar meson octet of pions, kaons and the eta meson. This mass gap is in turn comparable 
to 471"/^, the scale associated with spontaneous chiral symmetry breaking in QCD. 



3 Chiral Effective Field Theory 

Low-energy QCD is the physics of systems of light quarks at energy and momentum scales smaller than 
the 1 GeV mass gap observed in the hadron spectrum. This 1 GeV scale set by 47rf n offers a natural 
separation between "light" and "heavy" (or, correspondingly, "fast" and "slow") degrees of freedom. 
The basic idea of an effective field theory is to introduce the active light particles as collective degrees 
of freedom, while the heavy particles are frozen and treated as (almost) static sources. The dynamics 
is described by an effective Lagrangian which incorporates all relevant symmetries of the underlying 
fundamental theory. In QCD, confinement and spontaneous chiral symmetry breaking implies that the 
"fast" degrees of freedom are the Nambu-Goldstone bosons. With Yukawa's pion in mind, we restrict 
ourselves to Nf — 2. 



3.1 The Nambu-Goldstone Boson Sector 

We first briefly summarize the steps [8] [9] required in the pure meson sector (baryon number B = 0) 
and later for the pion-nucleon sector (B = 1). A chiral field is introduced as 

U(x) = exp (1. r ■ n(x)) G SU(2) , (8) 

with the Goldstone pion fields tt(x) divided by the pion decay constant f n in the chiral limit (m n — >■ 0). 
The QCD Lagrangian is replaced by an effective Lagrangian which involves the chiral field U(x) and 
its derivatives: 

£ Q CD^£eff(U,d»U,..:). (9) 

Goldstone bosons interact only when they carry non-zero four-momenta, so the low-energy expansion 
of OH]) is an ordering in powers of d^U. Lorentz invariance permits only even numbers of derivatives. 
One writes a series 

£ eff = 4 2 i + 4 4 2 + -, (io) 

where the leading term (the non-linear sigma model) involves just two derivatives: 

£g = S-triVUdrtf) . (11) 

At fourth order, the terms permitted by symmetries are (apart from further contributions involving the 
light quark mass m q and external fields, not shown here): 

£($ = l ± kr^Ud^)} 2 + ^(dJJdyrf) ti{d»Ud v U ] ) + ... , (12) 
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The constants £±,£2 (following canonical notations [9J) absorb divergences of loops, and their finite 
scale-dependent parts must be determined from experiment. 

The symmetry breaking mass term is small, so that it can be handled perturbatively, together with 
the power series in small momenta. The leading contribution introduces a term linear in the quark 
mass matrix m q = di&g(m u ,m d ): 

4 2 i = Str^C/cWt) + ll m l tr(f/ + [/t) ; (13) 

with m 2 ~ (m u + ma). The fourth order term 41r a ^ so receives explicit chiral symmetry breaking 
contributions (proportional to m q and m q ) with additional low-energy constants £3,4. 

To the extent that the effective Lagrangian includes all terms allowed by symmetries of QCD, 
the chiral effective field theory is the low-energy equivalent J2TJ [22] of the original QCD Lagrangian. 
Given the effective Lagrangian, the framework for systematic perturbative calculations of (on-shell) 
S- matrix elements involving Goldstone bosons, Chiral Perturbation Theory (ChPT), is then defined by 
the following rules: 

Collect all Feynman diagrams generated by C e ff- Classify all terms according to powers of a small 
quantity p which stands generically for three-momenta or energies of Goldstone bosons, or for the pion 
mass m n . The small expansion parameter is q/(47rf n ). Loops are evaluated in dimensional regular ization 
and get renormalized by appropriate chiral counter terms. 
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3.2 The Chiral Pion-Nucleon Effective Lagrangian 

The prominent role played by the pion as a Goldstone boson of spontaneously broken chiral symmetry 
has its impact on the low-energy structure and dynamics of nucleons as well [23] ■ When probing the 
nucleon with long- wavelength electroweak fields, a substantial part of the response comes from the pion 
cloud, the "soft" surface of the nucleon. The calculational framework for this, baryon chiral perturbation 
theory [2U [25] has been applied quite successfully to a variety of low-energy processes (such as low- 
energy pion-nucleon scattering, threshold pion photo- and electroproduction and Compton scattering 
on the nucleon). 

Consider now the sector with baryon number B = 1 and the physics of the pion-nucleon system. 

The nucleon is represented by an isospin-1/2 doublet, Dirac spinor \l/ = ( j of proton and neutron. 

The free Lagrangian 

4? e = *(z 7/ ^ - M )* (14) 

includes the nucleon mass in the chiral limit, Mq. One should note that the nucleon, unlike the pion, 
has a large mass of the same order as the chiral symmetry breaking scale 47r/ ir , which survives in the 
limit of vanishing quark masses, m u> d — > 0. 

The previous pure meson Lagrangian C e ff is now replaced by C e ff(U, d^U, *, ...) which also includes 
the nucleon field. The additional term involving the nucleon, denoted by C^n, is expanded again in 
powers of derivatives (external momenta) of the Goldstone boson field and of the quark masses: 

C wN = £ 7rJV + C nN + ... (15) 

In the leading term, C %N there is a replacement of <9 M by a chiral covariant derivative which introduces a 
vector current coupling between the pions and the nucleon. Secondly, there is an axial vector coupling. 
This structure of the irN effective Lagrangian is again dictated by chiral symmetry. We have 

C%1 = * [ii^ + r") - M + g A l^ u»] * , (16) 

with vector and axial vector quantities involving the Goldstone boson (pion) fields via £ = \/U in the 
form: 

^ = ^ t ,^] = i ^r-(vfx^7?) + ... , (17) 

«" = ji? ,&*£} = —?•&>*+... , (18) 

where the last steps result when expanding T^ and u^ to leading order in the pion fields. So far, the 
only parameters that enter are the nucleon mass M , the pion decay constant /„., and the nucleon axial 
vector coupling constant g^, all three taken in the chiral limit. 

At next-to-leading order, C n ^, the chiral symmetry breaking quark mass term enters. It has the 
effect of shifting the nucleon mass from its value in the chiral limit to the physical one. The nucleon 
sigma term 

a N = m JL = (NlmJuu + dd)\N) (19) 

om q 

measures the contribution of the non- vanishing quark mass, m q = \{rn u -\-mf), to the nucleon mass M N . 
Its empirical value is in the range a^ — (45 ± 8) MeV and has been deduced [26J by a sophisticated 
extrapolation of low-energy pion-nucleon data using dispersion relation techniques. Up to this point, 
the ttN effective Lagrangian, expanded to second order in the pion field, has the form 

■"N _ vfr('„-_, nV 1\/T Wr ^ A <Tr. 

2/. 



C N eff = *(i 7A4 d" - Mv)* - |ft 7/ ,7 5 r* • 3"7? 
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Figure 3: Hierarchical organization of nuclear forces in chiral effective field theory, 
where we have not shown a series of additional terms involving (<9 M 7r) 2 arising from the complete 

(2) 

Lagrangian C^jy. These terms come with further low-energy constants 03,4 encoding physics at smaller 
distances or higher energies. These constants need to be fitted to experimental data, e.g. from pion- 
nucleon scattering. 

The "effectiveness" of such an effective field theory relies on the proper identification of the active 
low-energy degrees of freedom. Pion-nucleon scattering is known to be dominated by the p-wave A (1232) 
resonance with spin and isospin 3/2. The excitation energy of this resonance, given by the mass 
difference A = Ma — Mn ~ 293 MeV, is not large, although it does not vanish in the chiral limit. If 
the physics of the A (1232) is absorbed in low-energy constants of an effective theory that works with 
pions and nucleons only (as done in heavy-baryon ChPT), the limits of applicability of such a theory is 
clearly narrowed down to an energy-momentum range small compared to A. The B = 1 chiral effective 
Lagrangian is therefore often extended [27] by incorporating the A(1232) isobar as an explicit degree 
of freedom. 



4 Chiral Nuclear Interactions 



4.1 Two-nucleon interaction 



In this section, we briefly review the structure of the chiral nucleon-nucleon potential in momentum- 
space [10J IH [12]. It amounts to a series of contributions from explicit multi-pion exchanges and 
short-distance contact-terms (shown in Fig. [3]) controlled by chiral symmetry and an expansion in small 
external momenta. The leading-order contribution consists of the well-known one-pion exchange piece 
and two contact-terms operating in the spin-singlet and spin-triplet S"-waves: 



V, 



(LO) 



9% 



ai-q a 2 -q 



NN 



4/2 m 2 +g 2 



f\ ■ r 2 + Cs + Ct Ox'B%. 



(21) 



Here, 0\p and 7^,2 denote usual spin- and isospin operators of the two nucleons and q is the momentum 
transfer between both nucleons. At next-to-leading order the two-pion exchange diagrams generated by 



the vertices of the chiral 7riV-Lagrangian C w ^ come into play. The corresponding potential reads 

^" LO) = 38i^j{ 4m - (1 + iA ~ iA) + qHl + 10A - 23A) ~ I^T?} L{q) ?1 ' ?2 

+ -^TfI L (l) fa • °* <? ~ *i ■ # ** • 9 ) > ( 22 ) 

with the logarithmic loop function 



,/4mJ + g 2 g + J4mj + g 2 

L(q) = l -— In- ¥ ^— . (23) 

W q 2m^ K J 

Notably, there are only contributions to the isovector central channel and the isoscalar spin-spin and 
tensor channel of the NN-interaction. Additional polynomial pieces generated by the pion-loops have 
been left out in the above expression. Without loss of information, these can be absorbed into the most 
general contact-term at NLO: 

+C 6 <?i ■ q a 2 ■ q + C 7 a 1 ■ k a 2 ■ k , (24) 

with k = (p + p')/2 the half-sum of initial and final nucleon momenta. The low-energy contants 

C*i, . . . , C-j are adjustable parameters to be determined in fits to empirical NN-phase shifts. Note that 

C*5 parametrizes the strength of the short-distance spin-orbit NN-interaction. In fact it dominates over 

the finite-range contributions to the spin-orbit interaction which arise at higher orders. The low-energy 

constant C5 will reappear in the discussion of the nuclear energy density functional, where the same 

coupling determines the (self-consistent) single-particle spin-orbit potential in finite nuclei proportional 

to the density gradient. 

The most important pieces from chiral two-pion exchange which generate attraction in the isoscalar 

central channel and reduce the too strong l7r-exchange isovector tensor force are still absent at this level. 

(2) 
These arise first from subleading 27r-exchange through the chiral ttttNN contact couplings 01,3,4 m ^n 

or from the inclusion of the A(1232)-isobar as an explicit degree of freedom. In the first approach one 

gets [28]: 



t/ (2tt-n 2 lo) 3g 

V NN 



2 



9 2 A c 

327T/4 



c 3 <f + 2m^(c 3 - 2ci) (2ml + q)Mq) 



2 
+Jr%( 4rr ^ + <l 2 ) A (<l) (^1 ■a 2 q 2 -a 1 -qa 2 -q)r 1 -r 2 , (25) 



with the loop function 



1 q 

A(q) = — arctan . (26) 

w 2q 2m n v ; 

At the same order there are additional relativistic 1/Mjv- correct ions to 27r-exchange. Their explicit 
form depends on the precise definition of the nucleon- nucleon potential Vaw, which by itself is not an 
observable. Covariant perturbation theory [28] and the method of unitary transformations [10] thus 
lead to slightly different expressions for these small l/M7v-corrections. As the state of the art, the 
chiral NN-potential has been constructed up to order N 3 LO and it includes two-loop 2-7r-exchange, 3n- 
exchange and contact-terms quartic in the momenta parametrized by 15 additional low-energy constants 
D\, . . . , D15. When inserted into the Lippmann-Schwinger equation (in order to solve for the unitary S- 
matrix or the T-matrix) the chiral NN-potential is multiplied by an exponential regulator function with 

9 



a cutoff scale A = 500 — 700 MeV in order to restrict the potential to the low-momentum region where 
chiral perturbation theory is applicable. Furthermore, methods of spectral function regularization [TU] 
have been employed in order to eliminate the high-momentum region in the pion-loop integrals directly. 
In this case the loop functions L(q) and A(q) receive an additional dependence on a regulator scale A. 
At order N 3 LO the chiral NN-potential reaches the quality of a "high-precision" potential in repro- 
ducing empirical NN-phase shifts and deuteron properties. At the same time it provides the appropriate 
two-body interaction constrained by chiral symmetry of QCD for nuclear few- and many-body calcula- 
tions. 

4.2 Nuclear Many- Body Forces 

Within the chiral effective field theory framework employing nucleons and pions as the explicit degrees 
of freedom, the leading-order contribution to the nuclear three-body force arises at order N 2 LO and 
consists of three terms. The two-pion exchange three-nucleon force contains terms proportional to the 
low-energy constants c±, C3, and C4 and has the form 

t/ (2tt) _ v- 9 A jj " gj &j " 9 j J?<*8a r (97 \ 

3N ~ifcwtf+ macs' +™}) F «* T,,i ' (7) 



where qt denotes the difference between the final and initial momenta of nucleon i, and the isospin 
tensor F°^ ]k is explicitly written 

Fg = 5 a ? (-4 Cl m 2 + 2c 3 £ • q 3 ) + c 4 e^r fe 7 a k • (gj x q 3 ) , (28) 



resulting in C\ and C3 terms proportional to f; ■ f ; and the C4 term proportional to T% ■ (7$ x Tj). The 
low-energy constants ci, C3, and C4 can be fit to pion- nucleon [29] or nucleon- nucleon scattering data 

[301 EH- 

The one-pion exchange term in the N 2 LO chiral three-nucleon interaction is proportional to the 
low-energy constant cjj and has the form 

while the three-nucleon contact interaction introduces the low-energy constant c^: 

V$= E o|a^^ (3°) 

where A = 700 MeV sets a natural scale. The low-energy constants cd and ce, associated with the 
one-pion exchange and contact three-nucleon forces respectively, ideally are fit to properties of three- 
body systems only. The binding energies of A = 3 nuclei and the /3-decay lifetime of 3 H provide largely 
uncorrected nuclear observables with which to fit these two low-energy constants |32j. 

At order N 3 LO in the chiral power counting, corresponding to 0(q A ) in powers of momentum, 
additional three- and four-nucleon force contributions arise without any additional undetermined low- 
energy constants. The N 3 LO three-body force can be written schematically as 

where the individual terms denote specific topologies. The 2tt and In — cont. topologies are present 
already in the N 2 LO chiral three-nucleon force, whereas the In — In, ring, and In — cont. topologies 
enter first at order N 3 LO. All contributions have been worked out and presented in refs. [331 |3U 135] . 
The N 3 LO four-nucleon force requires the evaluation of much fewer diagrams than the N 3 LO three-body 
force, and the resulting analytic expressions (presented in ref. [36]) are also considerably simpler. 
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Figure 4: Left: Total n + p cross section in the region of the A(1232) resonance (adapted from |37j). 
Right: Difference of polarized Compton scattering cross sections of the proton, 03/2 and 01/2, referring 
to channels with total angular momentum | and |, respectively. Data taken from |38j . Curves represent 
dispersion relation and multipole analysis cited in [38] . 

4.3 Role of Explicit A (1232) Degrees of Freedom 

The standard version of chiral meson-baryon effective field theory works with pions and nucleons only, 
both of which are stable particles with respect to the strong interaction. Effects of the A(1232) are 
encoded in low-energy constants such as C3 and C4 that are determined either by pion-nucleon scattering 
data or by fits to nucleon-nucleon phase shifts. 

On the other hand, the mass difference between the A(1232) isobar and the nucleon is only about 
0.3 GeV ~ 2m n , yet another "small" scale compared to the chiral symmetry breaking scale, 4^/^ ~ 
1 GeV. This suggests incorporating the A(1232) as an additional explicit degree of freedom in the 
effective field theory. In fact the A isobar is by far the dominant feature in the excitation spectrum of 
the nucleon observed in pion and photon scattering measurements. An example is the strong spin-isospin 
excitation observed in pion-nucleon scattering (see Fig. HJ left). The A(1232) is also seen prominently 
in polarized Compton scattering on the proton (Fig. HI right). 

In the low-energy expansion of the spin-independent nN scattering amplitude, the strong spin-isospin 
response of the nucleon manifests itself in a large "axial" polarizabilty: 



a 



(A) 



9\ 



M, 



5fm 3 



(32) 



N , 



Here the factor g\/f% comes from the axial vector coupling of the pion that drives the ttN — » A 
transition. The mere magnitude of this polarizability, several times the volume of the nucleon itself, 
already illustrates the importance of this effect. 

When implemented in the nucleon-nucleon interaction, the A isobar plays an important role in two- 
pion exchange processes such as the one shown in Fig. [5] (left). This mechanism contributes a large 
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Figure 5: Left panel: two-pion exchange involving a virtual N — y A — y N transition. Right panel: 
three-nucleon interaction generated by the same mechanism. 



part of the attractive isoscalar central part of the NN interaction [39J, the one often parametrized in 
terms of an ad-hoc "sigma" boson in phenomenological one-boson exchange potentials. A parameter- 
free calculation of this isoscalar central potential with single and double A excitation [40] agrees almost 
perfectly with phenomenological "a" exchange at distances r > 2 fm, but not at shorter distances. The 
detailed behavior of the 27r-exchange isoscalar central potential with virtual excitation of a single A is 
instead more reminiscent of a van der Waals potential: 



Vi 



NAi 
C \ 



3g\a 



(A) 



-2m- K r 



(8tt/ ot 



-P(m^r) 



(33) 



where P(x) = 6 + 12x + 10x 2 + 4x 3 + x 4 is a fourth-order polynomial in x 



m n r. The familiar r~ 6 
dependence of the non-relativistic van der Waals interaction emerges in the chiral limit, m n —y 0. 

When "integrating out" the A degrees of freedom, the van der Waals - like two-pion exchange 
mechanism just described - gives rise to an important effective three-nucleon interaction, Fig. 
ref2piDelta (right), that was actually suggested already more than half a century ago by Fujita and 
Miyazawa [H]. In a many-body framework with explicit A degrees of freedom, such a mechanism 
automatically occurs as an iterated one-pion exchange process involving the N — y A — y N transition. 
In this case the low-energy constants c 3 and c 4 in eq. ff28|) . the ones related to p-wave pion-nucleon 
scattering, are readjusted and considerably reduced in magnitude since they have to account only for 
the remaining non-resonant background, once the A(1232) is treated explicitly. Then important physics 
of the A such as the terms illustrated in Fig. (left) are actually promoted from N 2 LO to NLO in the 
chiral hierarchy of NN interaction terms [TD] , leading to improved convergence. 

These considerations about the role of the A(1232) in the low-energy description of nuclear forces 
should already suffice to answer one of the principal questions at the starting point of the chiral effective 
field theory approach: how much information about the intrinsic structure of the nucleon is actually 
needed in order for this theory to work efficiently? The A(1232) is indeed by far the dominant feature 
of nucleon structure in the relevant low-energy range. Other, less prominent properties such as baryon 
resonances appearing at higher energies are conveniently absorbed in the set of low-energy constants. 



5 Nuclear Chiral Dynamics 

5.1 In-medium Chiral Perturbation Theory: Nuclear Matter 

The tool to investigate the implications of spontaneous and explicit chiral symmetry breaking in QCD 
is chiral perturbation theory [9]. Observables are calculated within the framework of an effective field 
theory of interacting pions (Goldstone bosons) and nucleons [25] - The diagrammatic expansion in the 
number of loops has a one-to-one correspondence to an expansion in small external momenta and the 
pion (or light quark) mass. In nuclear matter, the relevant momentum scale is the Fermi momentum 
kf, related to the nucleon density by p = 2k 3 J3ir 2 . At the empirical saturation density po — 0.16 fm -3 
the Fermi momentum and the pion mass are of comparable magnitude, fe/o — 2m,,-. This immediately 
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Figure 6: One-pion exchange Fock diagram and iterated one-pion exchange Hartree and Fock diagrams 
contributing to the energy density pE(kf) of isospin-symmetric nuclear matter. 

implies that pions must be included as explicit degrees of freedom: their propagation in nuclear matter 
is resolvable and thus relevant at the densities of interest. With two small scales, kf and m n , at hand 
the nuclear matter equation of state as obtained from in-medium chiral perturbation theory will be 
given by an expansion in powers of the Fermi momentum. The expansion coefficients are non-trivial 
functions of kf/m n , the dimensionless ratio of the two relevant small scales inherent to the problem. 

The only new feature in performing calculations in a nuclear many-body system (as compared to 
scattering processes in the vacuum) is the in-medium nucleon propagator. For a non-relativistic nucleon 
with four-momentum p M = (po,p) it reads: 



G(po,p] 



p - p 2 /2M N + ie 



2w8(p -p 2 /2M N )9(k f -\p\) 



(34) 



where the second term, the so-called "medium insertion" , accounts for the presence of a filled Fermi 
sea of nucleons. The expression for G(po,p) is identical to the conventional decomposition of the in- 
medium propagator into particle and hole components used in non-relativistic many-body perturbation 
theory. With the decomposition into vacuum propagator and medium insertion, closed multi-loop 
diagrams representing the energy density can then be organized systematically in the number of medium 
insertions. The non-relativistic approximation is also compatible with the fcj-expansion at leading order. 
We outline the leading contributions to the energy per particle E{kf) of isospin-symmetric spin- 
saturated nuclear matter [42] . The relativistically improved kinetic energy is of course well-known: 



E(k f ) 



(kin) 



10M, 



N 



3k f 

56M^ 



(35) 



Nuclear chiral dynamics up to three loop-order introduces the closed one- and two-ring diagrams shown 
in Fig. O The contribution of the one-pion exchange Fock term (represented by the left diagram) to 
the energy per particle reads: 



E(k f ) 



(i ff ) 



(47T/,) 



m 3 1 3u / 3 

1 h arctan 2u — [ h 

3 8u 4 ^Su 



32m 3 



ln(l + 4m 2 



(36) 



with the dimensionless variable u = kf/m n . Next in the /^-expansion come the contributions from 
iterated (second-order) l7r-exchange. The latter refers to the (two-particle) reducible part of the planar 
27r-exchange box diagram which gets enhanced by a small energy denominator proportional to differences 
of nucleon kinetic energies. The corresponding Hartree term (see middle diagram in Fig. [6]) with two 
medium insertions reads: 



E{k f ) 



{m) = 39 1 m ^7: \t- - 59u + ( 6o + 32 " 2 ) arctan2 « - ( A + Pi in (! + ^ 

b(8Tt) 6 j£ V2u v 8w d 2u> 



(37) 



The analytical evaluation of this contribution is based on the iterated l7r-exchange amplitude in the 
forward direction. For the Fock diagram (see right diagram in Fig. [6]) with two medium insertions 
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one encounters, in the same way, the iterated l7r-exchange amplitude in the backward direction. The 
corresponding contribution to the energy per particle reads: 



E(k f y 



{F2) _ g\M N mj / u 3 | /■« 3x(u-x) 2 (2u + x) 



Wf* 



— + / dx 



2w 3 (l + 2x 2 ) 



1 + 8x + 8a: ) arctanx — (1 + Ax ) arctan2x 



(38) 

Note that the expressions for these two-body terms do not include the contribution of a linear divergence 
J °° dl 1 of the momentum-space loop integral. Such a linear divergence is set to zero in dimensional 
regularization, but employing a cutoff regularization gives a contribution proportional to the momentum 
cutoff A. There is an interpretational problem with results of dimensional regularization. Contributions 
which are expected to be attractive according to second-order perturbation theory arguments show only 
their finite repulsive parts. We restore the attractive component via a term linear in the cutoff and the 
density 

10g\M, 



E{k f ) 



(A) 



<JV 



A h 3 



(39) 



to which the Hartree and Fock diagrams have contributed in the ratio 4:1. 

The Pauli blocking corrections due to the nuclear medium are included through diagrams with three 
medium insertions. The corresponding Hartree term has the form: 



E(k f ) 



(H3) 



9g\M N ml 
(4vr/ 7r )%3 J 



KAJiXj iXj 



dy 



2uxy + (u — x y ) In 



u + xy 
u — xy 



2s 2 + s 4 
1 + s 2 



21n(l + s 2 ) , (40) 



with the abbreviation s = xy + y/W 
medium insertions has the form: 



x 2 y 2 . Similarly, the Fock-diagram in Fig. |6] containing three 



E{k f ) 



(F3) 



9g\M N mi 



(AttJ^u 3 Jo 
with the auxiliary function 



dx 



G 2 



+ 



x 



dy 



1 dz yz8{y 2 + z 2 - 
-l \yz\yjy 2 + z 2 



s 2 -ln(l + s 2 ) ln(l + t 2 )-t 



G 



u{l + u 2 + x 2 )-hl + {u + x) 2 ][l + {u-x) 2 }\n \ + { ' l + X) l 
4x 1 + (u - x) 2 



(41) 
(42) 



and t = xz + y/u 2 — x 2 + x 2 z 2 . Let us emphasize that Pauli blocking has been treated here exactly, 
without employing any simplification due to angular averaging. 

Summing the contributions to the energy per particle E{kf) listed above and treating the cutoff scale 
A as an adjustable parameter to incorporate (in a global fashion) unresolved short-distance dynamics, 
one obtains the equation of state of nuclear matter as shown in Fig. [71 With a fine-tuned cutoff 
scale of A = 611 MeV, the saturation minimum lies at po — 0.173 fm~ 3 and E(kf ) = — 15.3MeV. 
The nuclear matter compressibility, related to the curvature of the saturation line at its minimum, 
comes out as K = k 2 E"(kf ) = 252 MeV, in good agreement with the presently accepted empirical 
value of K = (250 ± 50) MeV [13, HU 03]. Although this looks at first sight to be a very successful 
reproduction of nuclear matter bulk properties, one should not hide the fact that strong cancellations 
between individually large terms of opposite sign are involved here. For example, the adjusted term 
linear in density amounts to —177.4 MeV at saturation density. The feature that strong cancellations 
are involved, which is also common to other approaches to nuclear matter, does not allow to make 
definite statements about the convergence of the ^/-expansion or to give reliable theoretical error bars. 
The appropriate framework to address such questions is a many-body calculation based on chiral low- 
momentum interactions whose long- and short- distance parts are determined together by empirical NN 
phase shifts and not separated any further. 
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Figure 7: Saturation curve of nuclear matter obtained from l7r-exchange and iterated l7r-exchange 
together with a fine-tuned short- distance term linear in density. The dashed line stems from the many- 
body calculation of ref. 



Let us now give an explanation for how saturation of nuclear matter is achieved in the framework 
of in-medium chiral perturbation theory. For that purpose it is instructive to consider the following 
simple parametrization of the energy per particle 



3k 2 k 3 k A 

Bl - k '^mk- a 4 +fi 4' (43) 

which includes an attractive fc^-term and a repulsive fc^-term. This two-parameter form has generically 
a saturation minimum if a, (3 > 0. Its striking feature is than once a = 5.27 and /3 = 12.22 are adjusted 
to the empirical saturation point po = 0.16 fm~ 3 , Eq = —16 MeV the compressibility K ~ 240 MeV 
comes out correctly. Moreover, such a parametrized curve for E{kf) follows the results of sophisticated 
many-body calculations [16] up to quite high densities p ~ 1 fm -3 . 

In the chiral limit m n = the leading interaction contributions calculated from lir- and iterated 

'A 



l7r-exchange turn into exactly such a two-parameter form with the coefficient /3 of the ki-term. given by 



P = ^ (if) ( 47r2 + m - 24 ln 2 ) = 13 ' 6 > ( 44 ) 

where g n N = 9a^n I fn = 13.2 is the strong pion-nucleon coupling constant. This number is quite 
close to (3 = 12.22 as extracted from a realistic nuclear matter equation of state. The mechanism for 
nuclear matter to saturate can be summarized roughly as follows: while pion-exchange at second order 
generates the necessary attraction, the Pauli-blocking effects due to the nuclear medium counteract this 
attraction in the form of a repulsive contribution with a stronger density dependence (a ki-term). 

Calculations of nuclear matter in this framework have been extended further by including the (ir- 
reducible) two-pion exchange contributions in the medium [5H] ■ A compact form of the corresponding 
Fock term is given in terms of a (subtracted) spectral function representation: 

- I f°° ( Ak 3 

E(k f Y 27TF) = — / dn Im(V c + 3W C + 2/i V T + 6fi 2 W T ) 3fik f - -*- 

07T d J2m-n I oil 
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Figure 8: Hartree and Fock three-body diagrams related to 27r-exchange with virtual A(1232)-isobar 
excitation. They represent long-range interactions between three nucleons in the Fermi sea. 



8kj 



2k f 



Varctan^ + ^(m ? + ^)ln(l + 5)} 



(45) 



where ImV^, ImH 7 ^, ImVy and IhiVFt are the imaginary parts (or spectral functions) of the isoscalar 
and isovector central and tensor NN-amplitudes, respectively. At this level also the 27r-exchange with 
excitation of virtual A(1232)-isobars comes into play. It is known to provide a dominant mechanism for 
intermediate range attraction in the isoscalar central channel [ID]. Furthermore, the Pauli-blocking ef- 
fects to these A(1232)-induced 27r-exchange mechanisms, as represented by diagrams with three medium 
insertions, are equivalent to contributions from genuine long-range three-body forces. The two-ring 
Hartree diagram in Fig. [S] leads to the following contribution to the energy per particle 



E{k f ) 



(AH3) 



4 6 
9 A m l 

A(2tt^; 



2 1 

-u 6 + u 2 - 3m 4 + 5u 3 arctan2u - -(1 + 9m 2 ) ln(l + Au 2 ) 



(46) 



with A = 293 MeV the delta-nucleon mass splitting. On the other hand, the total contribution of the 
one-ring Fock diagrams in Fig. [8] reads: 



E{k f ) 



(AF3) _ 3g A m n u 



6„,-3 



dx 



4A(4tt/ u ) 4 Jo 
where we have introduced the two auxiliary functions: 

Aux 
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x- 
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1 + (u + x)' 
1 + (it -xf 



(49) 



Evidently, the three-body Hartree term is repulsive while the corresponding Fock term is weakly at- 
tractive. 

The subtraction constants asscociated with the spectral function representation eq. (T4"5"j) give rise to 
short-distance contributions of the form 



E{k f ) 



k" 



k<] 



k ( ) 



B 3 ^L + Bs-t + B 6 ' 



'M% 



'M% 



M% 



(50) 



where the last term stems from a three-body contact term. The parameters Bs,B^, Bq are assumed to 
represent all unresolved short-distance NN-dynamics relevant for isospin-symmetric nuclear matter at 
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low and moderate densities. They are adjusted to few empirical bulk properties of nuclear matter, and 
thereafter the approach is able to make predictions for more detailed quantities like the single-particle 
potential or the thermodynamic properties of nuclear matter at finite temperatures. Another benefit 
of this approach is that the pion-mass dependence of all interaction terms is explicitly known. By 
differentiating the nuclear matter equation of state with respect to the light quark mass or pion mass 
one gets access to the density and temperature dependent chiral condensate (qq)(p,T). 

5.2 Resummation Strategies for Neutron Matter at Low Densities 

Dilute degenerate many-fermion systems with large scattering lengths have attracted much interest in 
recent years [IT] . Experimental advances have allowed the possibility to tune (magnetically) atomic in- 
teractions in ultracold fermionic gases through so-called Feshbach resonances, enabling the exploration 
of weakly and strongly interacting many-body systems together with transitions from the supercon- 
ducting BCS phase to the Bose-Einstein condensed state. Of special interest is the so-called unitary 
limit in which the two-body interaction is capable of supporting a single bound state at zero energy and 
the S'-wave scattering length diverges, a — > oo. Under these conditions the strongly- interacting many- 
fermion system is scale invariant. The ground-state energy per particle is then determined by a universal 
number, the so-called Bertsch parameter £, which is equal to the ratio of the energy per particle of the 
strongly interacting system at unitary, E(kfp°°\ to that of a free Fermi gas, E(kfp ^ = 3k"j/10M. 

Low-density (p n = k 3 r/3ir 2 < 0.05 fm 3 ) neutron matter is expected to behave as a unitary Fermi 
gas due to the very large neutron- neutron S- wave scattering length, a nn ~ 19 fm [IS]. The results of a 
variety of sophisticated many-body calculations [19~| [50] indicate that at low densities the interactions 
behind the large neutron-neutron S"-wave scattering length lead to a reduction of the free Fermi gas 
energy by approximately a factor £ n „ ~ 1/2. 

Here, we outline the basic steps of a novel resummation technique [51] which allows to sum up the 
complete series of fermionic in-medium ladder diagrams generated by a contact interaction proportional 
to the scattering length a. The solution of this (restricted) problem is enabled through a different 
organization of the many-body calculation. Instead of treating particles and holes separately, one 
keeps them together and takes into account effects of the filled Fermi sea by a "medium-insertion". 
This approach is realized by the following identical rewriting of the (non-relativistic) in-medium single- 
particle propagator: 

G(P0 , P -) = ,( wi-y, + «*'-w 



p -p 2 /2M + ie p -p 2 /2M -ie, 

- 27r5(p -p 2 /2M) 9(k f - \p\) , (51) 



Po - p 2 /2M + ie 



for an internal fermion line with energy po and momentum p. In the above equation, M denotes the 
large fermion mass. 

Many-body contributions to the energy per particle come from diagrams with at least two medium- 
insertions. For the closed ladder diagram, shown to the left in Fig. [HJ this minimal pair of medium- 
insertions must be placed on opposing positions of the double-ring. By opening them one obtains a 
planar ladder diagram containing energy denominators that are all given as differences of fermion kinetic 
energies. It is only for this planar topology that the factors of M from the energy denominators balance 
the 1/M factors coming from the interaction vertices in such a way that in the non-relativistic limit a 
finite result remains (to any order a n ). The open ladder diagram, shown to the right in Fig. [HI comes 
in with additional medium-insertions on internal lines in all possible ways. Due to the special nature 
of the momentum-independent contact interaction, all multi-loop diagrams factorize (successive loops 
are independent), which allows their sum to be written in the form of a power of the in-medium loop. 
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Figure 9: Left: Closed multi-loop diagram representing a contribution to the energy density. Right: 
Planar ladder diagram obtained by opening at a minimal pair of adjacent medium- insertions. Wavy 
lines symbolize the contact interaction proportional to the scattering length a. 




The in-medium loop is decomposed into contributions from zero, one and two medium insertions: 
Bq + B\ + B 2 (see Fig. [TO]) . The first term B is the well-known rescattering in vacuum: 



B = Ana 
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(52) 



where the rule / °° dl 1 = of dimensional regularization has been exploited. In a regularization scheme 
employing a momentum-space cutoff, the emerging scale-dependent constant —2A/ti can be absorbed 
into a -1 , from which one defines the renormalized (or physical) scattering length. The expression for 
the contribution B\ to the in-medium loop with one medium insertion reads: 

Solving this integral, one obtains the following result for the real part: 

ak f 



(53) 
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(54) 



(55) 



S — K\ 'IS' 1 — S + K 

is a logarithmic function written in terms of two dimensionless variables s = \pi +P2I/2&/ and K = 
\pi — p 2 \/2kf. An additional constraint, s 2 + n 2 < 1, arises from the fact that both external momenta 
Pi and P2 lie inside the Fermi sphere. 

Finally, there is the contribution B 2 with two medium insertions which in fact is purely imaginary. 
Summing up all contributions gives for the total imaginary part of the in-medium loop: 

1111(50 + #i + B 2 ) = 4vra / -^- vr 5(l 2 - f) 



J (2-kY- 

i-e{k f -\p-i\)}[i-9{kf-\P+i\) 

+e(k f -\p-r\)e(k f -\p+r\)), 



(56) 
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Figure 10: In-medium loop with different numbers of medium insertions. The middle diagram has a 
reflected partner. The external momenta \pip\ < kf are inside the Fermi sphere. 

where terms with zero, one, and two step-functions have suitably arranged. The first term has the 
form [1 — #(...)] [1 — #(...)] and gives no contribution to the imaginary part since energy conservation 
prohibits the on-shell scattering of two particles within the Fermi sea into the region outside the Fermi 
sea. Consequently, the result for the imaginary part of the in-medium loop takes the form: 



lm{B + B 1 + B 2 ) = ^ = ak f I{s, k) 



(57) 



where I(s, k) is the non-smooth function 



I(s, k) = k 9(1 - s - k) + — (1 - s 2 - k 2 ) 6(k + s-1). 



(58) 



In this expression k is restricted to lie in the interval < k < a/1 — s 2 . Putting together the real and 
imaginary parts, the complex-valued in-medium loop is written 



Bn 



akf 



B x + B 2 = f -{R(s,K)-iirI{s,K)}. 



(59) 



We observe that if the contribution from the diagram with two medium-insertions is separated off, then 
the imaginary part of the expression in eq. ( )59l) changes sign: 



(60) 



ok 



B + B ± = f - {R(s, k) + i-n I(s, k)\ . 



Indeed, this particular property of the in-medium loop will play a critical role in the derivation of the 
correct expression for the energy per particle. 

In ref. [51] it was shown how the contributions to the energy per particle E(kf) at a given order in a n 
can be constructed from the in-medium loop and how the resulting series in akf can even be summed to 
all orders. A naive iteration method would suggest (R — inI) n ~ 1 as the integrand for the energy density 
at order a n , but it is complex- valued. The deficit of the naive iteration method lies in the fact that a 
closed diagram with a repeated double medium-insertion has a symmetry factor which is not included 
in the binomial series expansion of [(R + ml) + (— 2ml)] n ~ 1 . A detailed combinatorial analysis [5Tj 
shows that the j-th power of —2ml coming from the diagrams with repeated double medium-insertions 
has to be re- weighted by a factor l/(j + 1). This crucial observation leads to the following formula 



n-l 

3=0 



-2iniy 



n 



j + 1 2m In 



{{R + tirI) n -(R-t7rI) n } 



(61) 



where the right hand side is now manifestly real for all n. In this form the whole subset of ladder 
diagrams can be summed to all orders, since the series Z^LJ— akf(R ± i7il)/ir} n /n is solved easily 
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Figure 11: The energy per particle of neutron matter as a function of the neutron density p n . The 
dash-dotted line is a reproduction of the results from the variational many-body calculation of ref. 



in terms of a (complex) logarithm. Putting all pieces together the final expression for the resummed 
energy per particle reads 



E(k f ) 



(lad) 



vrM 



i VT 

ds s 2 



o 







<1kk arctan 



akr I[ 



S, K 



1 + 7r l akf R(s, k) 



(62) 



where the arctangent function refers to the usual branch with odd parity taking on values in the interval 

[-7T/2,7T/2]. 

The expansion of the resummed energy per particle E(kf)( lad > in powers of akf has been checked 
against known results from many-body perturbation theory based on the traditional particle-hole count- 
ing scheme J52j[53]. Up to the order a 4 where such results are presently available the agreement is perfect. 
Moreover, since only double integrals are involved, the pertinent coefficients can be calculated with very 
high numerical precision [5T] . 

As an application of eq. ( 16"2"j) for the exact resummation of in-medium ladder diagrams, one can 
consider the equation of state of pure neutron matter. The very large neutron-neutron scattering 
length a nn = (18.95 ± 0.40) fm |48j necessitates a non-perturbative treatment of neutron matter at low 
densities. In Fig. [11] we reproduce the energy per particle of neutron matter as a function of the density 
in the approximation E(kf)( laxi ) + 3fc?/10M [5TJ. The solid line is obtained from eq. (16"2"j) by substituting 
Q — O'nn — 18.95 fm for the scattering length and M = M n = 939.57 MeV for the large fermion mass. 
The dash-dotted line comes from a sophisticated variational many-body calculation [54] representative 
of realistic neutron matter calculations. In Fig.[TT]the dashed line corresponds to Steele's suggestion [53] 
of a simple geometrical series, E(kfp^ = — ak^[3M(n + 2akf)]~ 1 . Up to rather high neutron densities, 
p n ~ 0.2 fm -3 , where the dimensionless parameter a nn kf ~ 34, there is rather good agreement among 
the different calculations. For neutron densities beyond that of saturated nuclear matter, repulsive 
effects from three-body forces (present in the variational calculation of ref. [54J) begin to play a more 
important role. Including the neutron-neutron S-wave effective range r nn = (2.75 ± 0.11) fm can also 
modify the equation of state, as shown in ref. [55J. 
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Figure 12: Diagrammatic contributions to the quasiparticle interaction in neutron matter generated 
from the two-pion exchange, one-pion exchange, and contact three-nucleon forces. The short double- 
line denotes the summation over the filled Fermi sea of nucleons. Reflected diagrams of (2), (3), (4), 
and (5) are not shown. 

5.3 In- Medium Effective Nucleon-Nucleon Interaction 

In applications of chiral effective field theory to light nuclei, the N 2 LO chiral three-nucleon force is 
routinely implemented exactly [101 [56]. For heavier systems, however, the inclusion of three-nucleon 
forces is computationally demanding and requires resources approximately an order of magnitude larger 
than that for two-nucleon forces alone. An alternative strategy has therefore been explored in a number 
of recent calculations [57J EB1 [SS EH ED E2] of finite nuclei as well as infinite nuclear and neutron matter 
in which three-nucleon forces are replaced by medium-dependent effective two-body interactions. The 
three-body force in second-quantized notation 



V 3 n = - 5] (123\V\456)alalala6a 5 a4, 

^" 123456 



(63) 



with antisymmetrized matrix elements (123|V|456), is normal ordered with respect to a convenient 
reference state \Q) which yields 



V, 



3iV 



-^(ijk\V2, N \ijk) + -J2{ijl\V 3N \ijA) :a[a 4 : +- J2 <*12|Vajv|i45> :a\ala b a A : 



+ 
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36 



ijU 



J1245 



J2 (123|V 3Ar |456) \a\alalaQa b a 4: 



(64) 



123456 



where the (alphabetic) indices i, j, k represent filled orbitals in the reference state, and : : denotes the 
normal-ordered product of operators satisfying 



:a\ . . . a n : \Q) = 0. 



(65) 



In refs. [581 EH] the reference state was chosen to be a filled Fermi sea of noninteracting nucleons 
characterized by the density as well as the isospin asymmetry parameter 5 np = (p n — p p )/(p n + Pp)- 
With the additional kinematical restriction that the scattering takes place in the rest frame of the 
medium, the on-shell scattering amplitude (equivalent to the normal-ordered two-body part of the full 
three-nucleon force) has the same form as the free-space nucleon-nucleon interaction: 



Vr 



mcd 
NN 



V c + Wc n-T2 + [Vs + W s n ■ r 2 ] <?i • <7 2 + [V T + W T r x ■ f 2 ] a x ■ qa 2 ■ q 
+ [V LS + W LS ft ■ r 2 ] i(ai + a 2 ) ■ {qxp) + [Vq + Wq f x ■ f 2 ] a ± -(qx p)a 2 -(gxp), (66) 

where the different scalar functions V and W depend on p and q, the initial relative momentum and 
the momentum transfer, respectively. The resulting in-medium two-nucleon interaction can then be 
implemented in the many-body method of choice, provided that one takes care of the different symmetry 
factors in the zero- and one-body normal-ordered components of a three-nucleon force compared to those 
from two-body interactions [60] . 
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Figure 13: Contributions to the on-shell S- and .D-wave matrix elements from the density-dependent 
nucleon-nucleon interaction derived in ref . [58] EH] . The solid curve denotes the matrix elements of a 
renormalization-group evolved NN potential and the numbered lines show the modification due to the 
six classes of diagrams shown in Fig. 



For the N 2 LO chiral three-nucleon force, there are six topologically distinct diagrams contributing to 
the in-medium nucleon-nucleon interaction shown in Fig. [12] The contributions labeled (1), (2), and (4) 
give corrections to one-pion exchange, while contributions (5) and (6) adjust the strength of the nucleon- 
nucleon contact interaction. On the other hand, the Pauli-blocked two-pion exchange diagram, labeled 
(3) in Fig. [12] gives rise to central, spin-spin, tensor, spin-orbit, and quadratic spin-orbit contributions. 
In Fig. [13] we reproduce the on-shell momentum-space matrix elements of the renormalization-group 
evolved Idaho N 3 LO chiral two-nucleon force together with the modifications arising from the six dif- 
ferent components of the density- dependent nucleon-nucleon interaction in selected relative S and D 
waves. The values of the relevant low-energy constants are C\ = — 0.76GeV _1 , c 3 = — 4.78GeV~ 1 , 
C4 = — 3.96GeV _1 , Co = —2.06, and ce = —0.625. The dominant effects arise from two-pion exchange 
dynamics, but in all partial waves the pion self-energy correction V^y ' and the pion-exchange vertex 
correction V^n ' approximately cancel. This leaves Pauli-blocking effects in two-pion exchange, rep- 
resented by VftZ ' , together with additional repulsion from the chiral three-body contact interaction, 
encoded in V™n ' , as the dominant effects from three-body forces. 

The normal-ordering approximation, in which the residual three-nucleon force in eq. (164"]) is neglected, 
has proven to be a useful approximation for calculations in medium-mass nuclei where exact ab-initio 
many-body methods incorporating three-body forces are computationally challenging. Within this ap- 
proximation three-nucleon forces were shown to provide the microscopic origin of the anomalously-long 
beta-decay lifetime of 14 C [58J, which subsequent no-core shell model calculations [63J have confirmed 
and clarified. In Fig. [14] we show the Gamow- Teller strengths from low-lying states in 14 C to the ground 
state of 14 N computed with density- dependent low-momentum chiral nucleon-nucleon interactions. In 
fact, only the ground state to ground state transition receives significant medium modifications that re- 
sult in a strong enhancement of the 14 C lifetime. Further work has employed normal-ordered two-body 
interactions to study three-nucleon force effects on the neutron drip-line in oxygen isotopes [61] [64] 
as well as on the shell structure of calcium isotopes [621 ES]- Benchmark calculations of the normal 
ordering approximation in closed shell nuclei have been carried out in the importance-truncated no-core 
shell model, and it was found that beyond the lightest nuclei, normal-ordered Hamiltonians provide an 
accurate substitute for full three-nucleon forces [67] . In Fig. [15] we reproduce from ref. [67] the expec- 
tation values of the leading-order chiral three-nucleon force at different levels of the normal-ordering 
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Figure 14: Gamow- Teller strengths for transitions from the low-energy even-parity states of 14 C to the 
ground state of 14 N as a function of the nuclear density. 



approximation. For both 16 and 40 Ca, the difference between the expectation of the normal-ordered 
2B approximation, consisting of the first three terms in eq. f l64l . and the expectation value of the exact 
three-body interaction is negligible. 

Density-dependent nucleon-nucleon interactions have also been used extensively in calculations of 
the equation of state of isospin-symmetric nuclear matter and neutron matter [60l [68j [691 [70]. This has 
allowed the study of three- nucleon forces beyond the Hartree-Fock approximation, but certain topologies 
are necessarily omitted in such an approximation. The inclusion of the full second-order contribution to 
the energy per particle from a three-nucleon force has been carried out only in the case of a three-body 
contact interaction [71 j. There it was found that the complete set of second-order three-body diagrams 
gives a contribution to the energy per particle that is roughly half that from the density-dependent 
nucleon-nucleon interaction. 
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Figure 15: Ground-state expectation values of the N 2 LO chiral three-nucleon force at different normal- 
ordering approximations. The bars are labeled according to the highest n-body (nB) contribution 
included in the calculation. Similarity renormalization group evolved interactions were employed for 
two different values of the flow parameter a. Figure reproduced from ref. 
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Figure 16: Contributions to the real part of the on-shell (u = q 2 /(2M^)) momentum- and density- 
dependent optical potential from chiral two- and three-nucleon forces. The solid and dashed-dotted 
lines are the first- and second-order contributions from the two-body potential, while the dashed line is 
the first-order contribution from the chiral three-nucleon force. The vertical dotted line indicates the 
Fermi momentum, and the dashed-double-dotted line is the second-order contribution without three- 
body forces. 



5.4 Nuclear Mean Field 

The average single-particle potential acting on a nucleon in a finite nucleus is a central concept in nuclear 
structure and reaction theory. For negative energy bound states, the nuclear mean field is associated 
with the shell model potential, while for positive energy scattering states it is identified with the optical 
model potential. The latter is complex and strongly absorptive, in contrast to the shell model potential 
which is a real-valued quantity. While phenomenological optical potentials have been used extensively 
to describe reactions on target nuclei close to the valley of stability, microscopic optical potentials have 
no adjustable parameters and therefore provide a reliable basis for extending to reactions on exotic, 
neutron-rich nuclei that will be studied at the next generation of radioactive beam facilities. Numerous 
microscopic many-body methods have been used to compute the single-particle potential (self-energy), 
including the Brueckner-Hartree-Fock [72J [73J d] , Dirac-Brueckner-Hartree-Fock [TSJITB], and Green's 
function [77J, EH] methods. Here we review only the recent calculations [791 E3 EH] that have been 
carried out in the framework of in-medium chiral perturbation theory for infinite isospin-symmetric 
nuclear matter. 

The starting point is the high-precision Idaho chiral two-nucleon potential with a cutoff of A = 
500 MeV supplemented with the leading-order chiral three-nucleon force whose low-energy constants 
( Cl = -0.81 GeV-\ c 3 = -3.2 GeV- 1 , c 4 = 5.4 GeV" 1 , c D = -0.2, and c E = -0.205) have been fit 



to nucleon-nucleon scattering phase shifts [3T] and the binding energy and lifetime of the triton 
Within microscopic many-body theory, the nuclear mean field is identified with the nucleon self-energy, 
T,(f,f',E) = U(r,f',E) + iW(r,r',E). For a homogeneous medium the self-energy can be written 
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Figure 17: Left panel: Three-body Hartree-Fock contributions to the real part of the nuclear optical 
potential. The individual single-particle potentials from the three components of the N 2 LO chiral 
three-nucleon force are shown separately and plotted as a function of the momentum at nuclear matter 
saturation density. Right panel: Density dependence of the real part of the momentum-dependent 
optical potential at second order in perturbation theory from chiral two- and three-body forces. 



as a function of the momentum, energy, and density (or Fermi momentum) E(g, uj\ kf). The first two 
perturbative contributions from two-nucleon forces are given by 



E ( 2 j>f(q,u;kf) = ^{qhssrft^NlqhtSSxtt^m, 



,(2) 



E ( 2 fi(q,u;kf) 
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123 
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123 



\(piP3S 1 S 3 t 1 t 3 \V2 N \qh2SS2tt 2 ) 

w + e 2 - e x - e 3 + irj 
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-n 1 n 2 n 3 ( y 27i) 3 5(p 1 +p 3 -q- h 2 ), 



-n 1 n 2 n 3 (27i) 3 8(h 1 + h 3 -q-p2), (67) 



where V^n is the antisymmetrized two-body potential, rii = 9{kf — \k%\) is the zero-temperature Fermi 
distribution, rii — 1 — n %i and the sum is over the momentum, spin, and isospin of the intermediate 
states. The first-order contribution from three-body forces has the form 



,(i) 



S 3iV y(g, w; kf) = ^2(qhih 2 ; ssis 2 ; ttit 2 \V 3N \qhih2] ssis 2 ; ttit 2 }niTi2, 



(61 



12 



where V 3 n is the antisymmetrized three-body interaction. As discussed in the previous section, the 
density-dependent two-nucleon interaction [5H [59] can be used in the last equation of (|67|) to study the 
effects of three-nucleon forces beyond the Hartree-Fock approximation. 

In Fig. [TS] we reproduce results from ref. [SD] for the real part of the on-shell self-energy (u> = 
q 2 /(2Mjv)) as a function of the density and momentum. The leading-order Hartree-Fock contribution 
from two- and three-body forces are given by the solid lines and dashed lines, respectively Both 
contributions are real and explicitly energy independent, but they have qualitatively different properties. 
The Hartree-Fock contribution from two-body forces is attractive up to the maximum momentum 
plotted in Fig. dH] and decreases in magnitude as the momentum increases. The N 2 LO chiral three- 
nucleon force, on the other hand, gives a repulsive contribution to the single-particle potential that 
increases nearly linearly with the background density but which varies only mildly with the momentum 
of the propagating nucleon. As shown in the left panel of Fig. [TTJ, the two-pion exchange component 
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Figure 18: The imaginary part of the nuclear optical potential from chiral two- and three-body forces at 
second order in many-body perturbation theory. The vertical dotted line indicates the Fermi momentum 
k f . 



of the N 2 LO chiral three-nucleon force provides about 80% of the total repulsion and is responsible for 
nearly all of the momentum dependence of the three-body Hartree-Fock mean field. 

A notable feature of the contact and l7r-exchange three-body forces is that they generate a momentum- 
independent and a nearly momentum-independent single-particle potential, respectively. This implies 
a strong correlation in the low-energy constants cjj and ce, where variations along the line 



ce = a ■ Co + const 



(69) 



give nearly equivalent descriptions of the mean field. The constant of proportionality has the value 
a ~ 0.21 ± 0.02 and is weakly dependent on the momentum and density. As noted in ref. [80], in the 
chiral limit the single-particle potential associated with the one-pion exchange three-body force has a 
particularly simple form which then yields the correlation coefficient a = Qa/^ — 0.3. 

Second-order contributions to the single-particle potential from two- and three-body forces are shown 
as the dashed-dotted lines in Fig. [16] The momentum dependence of the second-order terms peaks close 
to the Fermi surface, and when combined with the Hartree-Fock mean field from two- and three-body 
forces the resulting momentum dependence is quite small below kf. The result is a nucleon effective 
mass at the Fermi surface that is close to the free-space mass. In Fig. [16] we show with the dashed- 
double-dotted line also the second-order contribution without three-nucleon forces. Second-order effects 
from three-body forces appear to be rather small and result in additional attraction at both low and 
high momenta. In the right panel of Fig. [T7] we show the total real part of the nucleon self-energy 
as a function of density and momentum. For low to moderate densities the single-particle potential is 
nearly momentum independent below kf. The overall well depth for scattering at zero incident energy 
(q = kf) is approximately — 57MeV at nuclear matter saturation density, which is within about 10% of 
the depth, — 52MeV, of phenomenological optical potentials [81J. 
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In Fig. [18] we show the imaginary part of the single-particle potential as a function of the density 
and momentum. At positive energies (measured with respect to the Fermi energy) the imaginary part 
of the potential is negative and strongly absorptive for particles propagating with large momenta. We 
see that the imaginary part vanishes quadratically close to the Fermi surface for both two- and three- 
nucleon force contributions, in agreement with Luttinger's theorem [82]. With two-body forces alone 
the imaginary part is approximately inversion-symmetric about the Fermi momentum, W(q, kf) ~ 
—W(2kf — q, kf), but three-nucleon forces lower the imaginary part at momenta well below and above 
kf. For intermediate-energy scattering, E ~ 100 MeV, the phenomenological strength of the imaginary 
part of the optical potential is approximately \W\ — 10 MeV [HI], which is much less than the value 
~ 30 MeV obtained in the second-order calculation with two- and three-body forces. Effective mass 
corrections [83, M], self-consistent single-particle energies in the denominators of eq. (loTj) . and an 
inclusion of an energy gap at the Fermi surface are all expected to reduce significantly the magnitude 
of the imaginary part of the calculated optical potential. Future comparisons to rare-isotope reactions 
will require also the inclusion of isospin asymmetry and the extension to finite nuclear systems. 



6 Density Functional Methods and Finite Nuclei 

In this section we review how chiral low-momentum interactions can be employed for the description 
of properties of finite nuclei [SSI ESI EH EBJ ESI ISS] • We leave aside the highly sophisticated ab-initio 
methods which have been developed over the past decades in order to solve (numerically) the nuclear 
A-body problem with a given two- and three-nucleon interaction. Here we concentrate on the nuclear 
energy density functional approach, which remains the many-body method of choice to compute the 
properties of medium- mass and heavy nuclei [9Tj 192] . Non-relativistic Skyrme functionals [93] 191] with 
a small number of adjustable parameters as well as relativistic mean-field models [951 ES] have been 
employed in such self-consistent mean-field calculations to describe a wide range of nuclear properties. 
An alternative and complementary scheme [SZ1 ESI ESj focuses less on fitting to experimental data and 
attempts instead to constrain the analytical form of the functional as well as the values of its couplings 
from a microscopic approach based on many-body perturbation theory with realistic two- and three- 
nucleon forces. The use of low-momentum nucleon-nucleon potentials [99j 1100] , rather than conventional 
hard-core interactions, is essential in this regard, because the former exhibit much improved conver- 
gence properties compared to the latter in perturbative solutions to the nuclear many-body problem. 
Indeed, when evolved low-momentum interactions are employed at second-order in perturbation theory, 
including also three-body forces, a good description of the bulk correlations in infinite nuclear matter 
|101[ [68] and in doubly-magic nuclei [102] can be achieved. 

Within the framework of many-body perturbation theory, contributions to the energy are given in 
terms of density-matrices convoluted with finite-range interaction kernels. The resulting expression 
for the energy is then highly non-local in both space and time. In order for mean-field calculations 
with such functionals to be numerically feasible in heavy open-shell nuclei, it is essential to develop 
approximations for these functionals in which only local densities and currents enter. For this purpose 
the density-matrix expansion is highly useful as it provides a means to remove the non-local character 
of the exchange contribution to the energy by expanding it in the form a generalized Skyrme functional 
with couplings that depend on the density. For many years such an approach has been based on the 
density- matrix expansion of Negele and Vautherin [103] , but recently Gebremariam, Duguet and Bogner 
[55] have suggested an improved version that specifically addresses the challenges that arise for spin- 
unsaturated nuclei. In fact, the phase-space averaging techniques they employed were shown to allow 
for a consistent expansion of both the spin-independent part as well as the spin-dependent part of the 
nuclear density-matrix. The improved features of the phase-space averaged density-matrix expansion 
have been studied [85J via the Fock energy densities of schematic finite-range central, spin-orbit, and 
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tensor interactions for a large group of semi-magic nuclei. 

Making use of these new developments a microscopically constrained nuclear energy density func- 
tional derived from the chiral NN potential at order N 2 LO has been presented in ref. [SB]- There it has 
been suggested that the density-dependent couplings associated with pion-exchange should be added 
to a standard Skyrme functional with re-adjusted parameters. In a subsequent study [15], it has been 
shown that this new energy density functional yields numerically stable results and that it exhibits 
a small yet systematic reduction of the \ 2 deviation in comparison to traditional Skyrme functionals 
without any explicit pion-exchange dynamics. 

In this section we report on a derivation of the nuclear energy density functional with improved 
(chiral) two- and three-nucleon interactions [89| [90] . For the two-body interaction the N 3 LO chiral 
NN-potential is used, which reaches at this order in the chiral expansion the quality of a high-precision 
nucleon-nucleon potential in reproducing empirical NN scattering phase shifts and properties of the 
deuteron. The N 3 LO chiral potential contains long-range one- and two-pion exchange interactions as well 
as a short-distance part that is parameterized in terms of 24 low-energy constants. The latter contact 
potential is written in momentum space and provides the most general contribution up to fourth power 
in nucleon momenta. This high-precision nucleon-nucleon potential is supplemented with the N 2 LO 
chiral three-nucleon force. In addition to the equation of state of infinite homogeneous nuclear matter, 
the energy density functional includes strength functions associated with the (Vp) 2 surface term and 
the spin-orbit coupling term. In phenomenological Skyrme parametrizations these strength functions 
are treated as constants, whereas in a microscopic approach based on realistic nuclear interactions the 
finite-range character of the pion-exchange terms gives rise to specific density dependences for these 
strength functions. 

6.1 Density- Matrix Expansion and Energy Density Functional 

The construction of an explicit nuclear energy density functional starts from the density-matrix given 
as a sum over the energy eigenfunctions \1/ Q (r ) associated with occupied orbitals of the non-relativistic 
many-body Fermi system. Gebremariam, Duguet and Bogner [85] have shown that it can be expanded 
in center-of-mass and relative coordinates, f and a, in the following way 



]T* a (r-a/2)*t(r + a/2) = j£. jl{akf) _ JL jl{akf) 



5 Pkf 4 



r - ^pkj - -V 2 p 



+ 2^rJi{ak f )a-{axJ) + ..., (70) 

where ji(x) = (sinx — xcosx)/x 2 is a spherical Bessel function. The right-hand side of eq. (ITU]) is 
written in terms of the nucleon density p(f) = 2/c 3 r (r*)/37r 2 = J2 a ^a(^*)^a(^*)> the kinetic energy 
density r{f) = S a V^(f) • V\l/ a (r ), and the spin-orbit density J{r) = i^ a W a {f)a x V x I / a (r ). 

Nuclear interactions derived within the framework of chiral effective field theory are generally given 
in momentum space. Therefore the Fourier transform of the expanded density-matrix eq. ( |70l) with 
respect to the coordinates a and r provides the appropriate tool for an efficient calculation of the 
nuclear energy density functional. This Fourier transform: 

F(p,q) = Jd 3 re- i ^! [ 9(k f -\p\) + ^[k f 5'(k f -\p\)-25(k f -\p\) 

x(r -\pk) -iv 2 p) ~^S(k f - \p\)ff ■ (px J)}, (71) 

generalizes the concept of "medium-insertion" to inhomogeneous many-nucleon systems characterized 
by the time-reversal-even fields p(r), r{r) and J{f). At a practical level T(p,q) extends the step- 
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function-like momentum distribution 0(kf — \p\) for infinite nuclear matter to inhomogeneous many- 
nucleon systems. Note that the delta-function 8(kf— \p\) in eq. (ITTj) gives weight to the nucleon-nucleon 
interactions in the vicinity of the local Fermi momentum \p\ = kf(r) only. 

Up to second order in spatial gradients, characterizing deviations from homogeneity, the energy 
density functional appropriate for N = Z even-even nuclei has the form 



£[p,T,J] = P E(p) + 



3 .„ir l fc 2 



in, " 



r - -pk f 



f 



+ Fr{ P ) 



2M N AM% 
+ (Vp) 2 F v (p) + Vp • JF so (p) + J 2 Fj(p) . (72) 

In this equation, E(p) represents the energy per particle of symmetric nuclear matter evaluated at the 
local nucleon density. The strength function F T (p) is associated with an effective density-dependent 
nucleon mass M*(p), and it is related to the nuclear single-particle potential U(p, kf): 

where p = 2kj-/37f 2 . The second equality gives the relation to the "spin and isospin independent" p- 
wave Fermi-liquid parameter f\(kf). A new feature of the improved density- matrix expansion is that it 
leads to the same concept of effective mass as established in Fermi-liquid theory for quasi-particles on 
the Fermi surface. For the original density-matrix expansion of Negele and Vautherin |103] this close 
relationship does not hold in general, e.g. when U(p, kf) deviates from a simple quadratic p-dependence. 
The strength function -Fv(p) associated with the (Vp) 2 surface term can be decomposed into two 
terms: 

Mp) = \^-+F d (p), (74) 

where Fd(p) is composed of all contributions for which the (Vp) 2 -factor originates directly from the 
momentum dependence of the interactions (expanded up to order q 2 and combined with a Fourier 
transformation). It is worth noting that only the nuclear matter component 9(kf — \p\) of the density- 
matrix expansion enters into the derivation of the strength function Fd(p). The next-to-last term in eq. 
( 1721) . namely Vp • J F so (p) with its associated strength function F so (p), gives the spin-orbit interaction 
in nuclei. The last term, J 2 Fj(p), in the expression for energy density is the quadratic spin-orbit term. 
It gives rise to an additional spin-orbit single-particle potential proportional to J, whereas the ordinary 
spin-orbit potential is provided by the density-gradient Vp. 

6.2 Two- and Three-Body Contributions 

Here, the two- and three-body contributions to the various density- dependent strength functions en- 
tering the nuclear energy density functional £[p, r, J] are presented. One would prefer to calculate the 
two-body contributions with perturbative nucleon-nucleon potentials, of which the low-momentum NN 
interaction VJ w-k [JOOJ is the prototypical example. Low-momentum two-body potentials are however 
non-local and given in terms of partial-wave matrix elements, which makes their implementation in the 
density-matrix expansion rather difficult. An explicit representation of the momentum-space nucleon- 
nucleon potential (without non-localities) in terms of spin- and isospin-operators is much more suitable 
for this purpose. As a convenient substitute for V[ow-k one can use the chiral NN-potential N 3 LOW 
developed in refs. [T2"j 1104] by imposing a sharp cutoff at the scale A = 414 MeV. This value of the 
cutoff coincides with the resolution scale below which evolved low-momentum NN potentials become 
nearly model-independent and exhibit desirable convergence properties in many-body perturbation the- 
ory calculations. The finite-range part of the N 3 LOW chiral NN-potential arises from one- and two-pion 
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exchange processes that have the general form: 



V, 



(*) 



NN 



Vc{q) + n-T 2 W c (q) + [V s (q) + n ■ r 2 W s (q)} 3\ • a 2 

+ [V T {q) + n ■ f 2 W T {q)] a x -qa 2 -q 

+ [Vso(q) + n-r 2 W so (q)] i{&i + ^2) • (q X p) , 



(75) 



where <f is the momentum transfer and pa single nucleon momentum. A special and convenient feature 
of V^n is that it is a local potential: all occurring functions Vc(q), • • • , Wso(o) depend only on the 
momentum transfer q, and quadratic spin-orbit components ~ <J\ ■ (q x p) a 2 ■ (q x p) are absent. 

In the Hartree-Fock approximation the finite-range part of the nucleon- nucleon potential arising from 
explicit pion exchange leads to the following two-body contributions to the energy density functional: 
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(76) 

(77) 

(78) 

(79) 

V T (q) + 3W T (q)]}, (80) 



where q = 2xkf. The above expressions are obtained from the density matrix-expansion employing the 
product of two medium-insertions r(pi, q ) T(p 2 , —q)). In eq. (!7H|) the double-prime denotes the second 
derivative, and it is worth noting that Fj(p) stays finite in the low-density limit kf — > 0. 

In addition to the finite-range parts of the NN potential there are two-body contributions from the 
zero-range contact forces. At order N 3 LO the corresponding expressions in momentum-space include 
constant, quadratic, and quartic terms in momenta |105j . The density-dependent strength functions 
arising from the contact NN potential read: 
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(82) 
(83) 
(84) 
(85) 



The 24 low-energy constants Cs,t, Gj and Dj are fit (with a sharp cutoff regulator at A = 414 MeV) 
to empirical nucleon-nucleon scattering phase shifts and properties of the deuteron |104] . It is worth 
mentioning that in the first-order Hartree-Fock approximation it is not necessary to include a regulator 
function since the nuclear interactions are probed only at small momenta |p lj2 | < kf < 285 MeV. The 
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Figure 19: Three-body diagrams related to the contact (eg) 
(ci, 02,03) components of the chiral three-nucleon interaction, 
medium insertion T(p, q) for inhomogeneous nuclear matter. 



l7r-exchange (cd), and 27r-exchange 
The short double-line symbolizes the 



regulator function as well as its associated cutoff scale A become relevant only at second order (and 
higher) . 

In the next step, the three-body contributions to the nuclear energy density functional S[p,r, J] 
are included. The leading order chiral three-nucleon interaction consists of a contact piece (with low- 
energy constant eg), a l7r-exchange part (proportional to c^), and a 27r-exchange component (with 
low-energy constants ci, c 3 and C4). In order simplify the treatment of the three-body correlations in 
inhomogeneous nuclear systems, it has been assumed [89J that the relevant product of density-matrices 
can be represented in the factorized form r(pi, q\) r(p2, Q2) r(p3, — q\ — qV) in momentum space. This 
factorization ansatz respects the correct infinite nuclear matter limit, but it involves approximations 
compared to more sophisticated treatments explored in ref. [88]. In fact, the present approach is similar 
to the method "DME-I" considered in ref. |88j . 

The three-body contact-term proportional to ce (see Fig. IT9~j) gives rise to a contribution to the 
energy per particle: 

c E kj 



E(p) 



6) 



12vr% 4 A ' 

that is quadratic in the density p = 2kjr/3n 2 . Obviously, all other strength functions F Tt d, S o,j(p) receive 
vanishing contributions due to the momentum-independence of the chiral contact interaction. The 
l7r-exchange component of the chiral 3N- interaction involves the low-energy constant cd- With three 
(inhomogeneous) medium insertions in the corresponding three-body diagram (shown in Fig. [T9l one 
finds the following expressions for the strength functions: 
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where u = kf/m^. We observe that there is no contribution to the spin-orbit strength F so (p), since 
neither the contact- vertex nor the l7r-exchange component produces a spin-orbit interaction. 

In the third diagram of Fig. [JJ5J we show the three-body contribution arising from the 27r-exchange 
Hartree term. With three (inhomogeneous) medium insertions, the corresponding energy density cou- 
pling strengths are written [89J: 



E(p) 



1 r 
9 A m l 

(27TUY 



(12ci — 10c 3 )m 3 arctan2w — -c 3 m 6 + 6(c 3 — Ci)u 4 



31 



+ (3ci - 2c 3 )u 2 + 



1 3m 2 

-(2c 3 -3 Cl ) + — (3c 3 -4 Cl ; 



ln(l+4w 2 



(91) 



Fr{ P ) 



1 4 
9A m % 

(2nf n y 



(5c 3 - Qci)u 2 + 



(c 3 - 2ci)w 2 



1 + Au 2 



+ 



2ci - -c 3 + 2(ci 



c 3 Jm 



ln(l + 4m 2 



(92) 



W 



2 

A TYb'fr I , 

-— — — -<^ (10ci - 23c 3 ) arctan2u + 16c 3 w 
(8tt) 2 /* 

, 7c 3 - 5ci 



» 



ln(l + 4m 2 ) + 



6c 3 u + 16(2c 3 — C\)v? 



3(1 + Au 2 y 



(93) 



Fso(p) 



Fj( P ) 



3gA m TT [2 



(87T) 2 /. 4 



H 



m„ 



57T 



|2f4 



■/7T 



-(4ci-3c 3 ) 
it 

3c 3 - 4ci 
it 



4c 3 m + 



It 



Cl) + 



2c 3 m + 



4m(2 Ci - c 3 ) 
1 + 4m 2 



3c 3 - 4ci 

2u 3 

4ci - 3c 3 

4m 3 



ln(l + 4m 2 



ln(l + 4m 2 



(94) 



(95) 



which depend only on the isoscalar coupling constants c\ and c 3 . Note that the expression for the 
spin-orbit strength in eq. (1941) provides the dominant part of the three-body contribution to F so (p), as 
suggested originally by Fujita and Miyazawa jH] in the context of A(1232)-resonance excitation. In 
the present chiral effective field theory treatment without explicit A resonances, the two-step process 
ttN —¥ A —¥ ttN is written as an equivalent nnNN contact vertex proportional to c 3 . 

Finally, we discuss the three-body 27r-exchange Fock contribution, which is shown as the fourth 
diagram of Fig. [HO This term consists of a single closed nucleon ring that generates for isospin- 
symmetric nuclear matter non-vanishing contributions from both the isoscalar and isovector inrNN 
contact vertices, the latter proportional to the low-energy constant C4. In contrast to the three-body 
27r-exchange Hartree diagram, the Fock contribution involves integrals (over three Fermi spheres) that 
cannot be computed in closed analytical form. The resulting rather lengthy expressions for the various 
strength functions can be found in section 4.4. of ref. [89J. 

Let us now discuss the results for the nuclear energy density functional obtained in the first-order 
Hartree- Fock approximation. In Fig. [20] we show the contributions to the energy per particle E(p) from 
chiral two- and three-body forces for densities up to p = 0.2 fm -3 . The dash-dotted line denotes the 
attractive two-body contributions, while the dashed line gives the repulsive three-body contributions. 
The results from the N 3 LOW chiral nucleon-nucleon interaction are compared to those from the universal 
low-momentum NN potential Vj ow -k obtained by integrating over its diagonal (on-shell) partial-wave 
matrix elements. Indeed, the treatment of the two-body interaction via the chiral potential N 3 LOW 
gives a fairly accurate reproduction of these results. The full line in Fig. [20] shows the sum of the two- 
and three-body contributions, which exhibits a first tendency for the saturation of nuclear matter. After 
including the kinetic energy E kin (p) = 3kj/10M N — 3kj/56Mff, however, the resulting saturation energy 
is still much too low. This particular feature of the Hartree- Fock approximation has been observed in 
similar studies |101[ 168] . A much improved description of the nuclear matter equation of state is achieved 
by treating the two-body (and three-body) interaction at least to second order. 

In Fig. EH we reproduce the contributions to the strength function F T (p) associated with the surface 
gradient term in the energy density functional. We note that the two-body parts derived from Vj ow -k 
and the chiral N 3 LOW potential are rather similar. The three-body part, which increases nearly linearly 
with the density, comes out relatively small. At nuclear matter saturation density po — 0.16 fm -3 it 
gives rise to a correction of about 20% compared to the two-body contribution. The bracketed term that 
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Figure 20: The energy per particle E(p) of isospin-symmetric nuclear matter derived from chiral nuclear 
interactions. 
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Figure 21: Contributions to the strength function F T (p) as a function of the nuclear density p. 

1 — , — , — , — , — i — , — , — , — , — i — , — , — , — , — i — , — , — , — , — i 



effective nucleon mass 




0.6 



0.05 0.1 0.15 0.2 

3, 



p[ftn ] 

Figure 22: Ratio of the effective nucleon mass M*(p) to the free nucleon mass M as function of p. 
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Figure 23: The strength function -Fy(p) of the surface term (Vp) 2 versus the nuclear density p. 
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Figure 24: Spin-orbit strength function F so (p) as a function of the density p. 
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Figure 25: Strength function Fj(p) of the squared spin-orbit density J 2 versus p. 
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multiplies the kinetic energy density r(f ) in the expression for the nuclear energy density functional in 
eq. (1721) can be interpreted as the reciprocal of a density-dependent effective nucleon mass: 



M*{p) = M N 



l-^ T +2M N F r (p)^ \ (96) 



In the improved density matrix expansion it is identical to the so-called "Landau" effective mass intro- 
duced in Fermi-liquid theory, and it is derived in the same way from the derivative of the single-particle 
potential U(p, kf) with respect to the momentum p at the Fermi surface p = kf. The relativistic increase 
in the mass is accounted for with the small correction term — H/2Mjy. In Fig. [22] we show the ratio 
of effective nucleon mass to the free nucleon mass M*(p)/Mjy as a function of the density p. In the 
Hartree-Fock approximation, the effective nucleon mass is reduced in comparison to the free-space mass, 
and the ratio reaches the value M*(p ) ~ 0.67M N at nuclear matter saturation density p — 0.16 fm~ 3 . 
This is consistent with the range 0.7 < M*(pq)/Mn < 1 determined from phenomenological Skyrme 
forces [931 Elj- On the other hand, second-order corrections from two-body forces significantly enhance 
the effective mass [106] . 

The strength function F\?(p) of the surface gradient term (Vp) 2 is plotted in Fig. [23] as a function 
of the nuclear density. At very low densities this coupling strength exhibits a pronounced increase in 
the two-body contribution (dash-dotted line) that has its origin in l7r-exchange. For larger densities, 
the two-body contributions are nearly independent of the density. The three-body contribution (dashed 
line) is sizable and negative. It reduces the two-body contribution such that the total value of -FV(p) 
(full line in Fig. 123]) decreases with increasing density p. To compare with phenomenology we also show 
the band (of constant F v (p)-values) spanned by parameterized Skyrme forces [93] |Mj. The Hartree- 
Fock result from realistic chiral two- and three-body forces is somewhat too small at densities close 
to po/2 = 0.08 fm~ 3 , where the surface energy in finite nuclei gains much of its weight. Iterated Itt- 
exchange has been considered in ref. |107] from which it may be deduced that a treatment of the 
two-body interaction to second order in perturbation theory will increase the values of Fy(p). 

Of special interest is the coupling strength F so (p) multiplying the spin-orbit term Vp ■ J. Contribu- 
tions from two- and three-body terms, as well as their total sum, are shown in Fig. [211 The coupling 
strength from two-body forces is dominated by the low-energy constant 3C 5 /8, as suggested also by the 
weak dependence of the dash-dotted line on the density p. Three-body forces induce a density-dependent 
spin-orbit interaction that considerably enhances the contribution from the N 3 LOW two-body force. 
The most significant contribution arises in the Hartree term (194]) proportional to the low-energy constant 
C3 = — 4.78GeV _1 . With this particular value of C3, the three-body spin-orbit strength is considerably 
larger than that proposed by Fujita and Miyazawa [UJ, where the A(1232)-excitation mechanism would 
correspond to c 3 = —g\/2A ~ — 2.9GeV _1 (A = 293 MeV is the delta-nucleon mass splitting). At 
densities around half that of saturated nuclear matter, po/2 = 0.08 fm -3 , the total Hartree-Fock con- 
tribution exceeds the empirical spin-orbit coupling strength F^ mp ^ (p) ~ 90 MeV fm 5 [931 El] by nearly 
50%. A compensating effect is therefore required, and it has been suggested in ref. [107] that the ln- 
exchange tensor force at second-order generates a spin-orbit coupling of the "wrong-sign" . It is then 
reasonable to assume that the low-momentum two-nucleon tensor potential, when treated to second- 
order in perturbation theory, will reduce the strength of the spin-orbit coupling F so (p) to a value close 
to that suggested by phenomenology. 

To conclude the description of the nuclear energy density functional in N = Z even-even nuclei, we 
show in Fig. [23] the strength function Fj(p) of the squared spin-orbit coupling. In contrast to all of the 
previous coupling strengths, Fj(p) receives only a very small contribution from three-body forces. One 
observes that the two-body contribution is, however, strongly density-dependent, and at low densities 
it reaches quite large values. This strong density dependence originates primarily from the l7r-exchange 
contribution, which for comparison we reproduce separately in Fig. [25] (dashed-double-dotted line). 
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At this point it should be emphasized that the J 2 -term in the energy density functional represents 
non-local Fock contributions from tensor forces, etc. It is therefore not surprising that there exists an 
outstanding l7r-exchange contribution to the strength function Fj(p). 

6.3 Isovector Part of the Nuclear Energy Density Functional 

The previous calculation of the nuclear energy density functional from chiral two- and three-nucleon 
forces can be extended in a straightforward manner to the isovector terms [90] relevant for nuclei with 
different proton and neutron densities. The additional isovector terms play an important role in the 
description of long isotopic chains of stable nuclei and for neutron-rich systems far from stability. By 
construction the density-matrix expansion of Gebremariam, Duguet and Bogner [HS] applies separately 
to proton and neutron orbitals. It is therefore easily adapted to the situation of isospin-asymmetric nu- 
clear systems with different proton and neutron local densities p P) „, r Pj „, J Pj „. After Fourier transforming 
to momentum space one obtains the medium insertion: 
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for the inhomogeneous isospin-asymmetric many-nucleon system. In the above expression, T3 denotes 
the third Pauli isospin matrix, and we have included only the relevant terms proportional to differences 
of proton and neutron densities: Pp — p n ,T p — T n , J p — J n . The two local Fermi momenta are related to the 
(particle) densities in the usual way: p p = k p /37i 2 , p n = kf l /3ir 2 , p — p p + p n — 2k 3 /3ir 2 . When working 
to quadratic order in deviations from isospin-symmetric nuclear systems, it is sufficient to employ an 
average Fermi momentum kf in the factors multiplying r p — r n and J p — J n . 

Up to second order in proton-neutron density differences and spatial gradients, the isovector part of 
the nuclear energy density functional reads: 

£iv[Pp, Pn, T p , T n , J p , J„] = ~(Pp- Pn) 2 A(p) + ~( T P ~ T n){Pp - Pn) G T {p) 

+iy Pp - Vp n ) 2 G v (p) + (V Pp - Vp n ) ■ (Jp - J n ) G so (p) + (J p - J n ) 2 Gj(p) . (98) 

Here, A(p) is the interaction part of the isospin- asymmetry energy of homogeneous nuclear matter. The 
non-interacting (or kinetic energy) contribution A^ n (p) = k"j/6M]y to the isospin-asymmetry energy is 
included in the energy density functional through the kinetic energy density term, £kin = (T p + T n )/2M^. 
The strength function G\?(p) multiplying the isovector surface term (Vp p — Vp n ) 2 can be decomposed 

as 

G v (p) = ^G T (p) + G d (p), (99) 

4p 

where Gd(p) denotes all those contributions for which the (Vp p — Vp„) 2 factor originates from the 
momentum dependence of the interactions in an expansion up to order q 2 . Performing a Fourier 
transformation converts this factor q 2 into (Vfc p — Vfc n ) 2 ~ (Vp p — Vp n ) 2 (ir/kf) 4: . The next-to-last 
term (Vp p — Vp„) • (J p — J n ) G so (p) in eq. (|98|) is the isovector spin-orbit interaction in nuclei. Depending 
on the sign and size of its strength function G so (p), the spin-orbit potentials for protons and neutrons 
receive different weightings from the gradients of the local proton and neutron densities. 
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Figure 26: Contributions to the asymmetry energy A(p) of nuclear matter. 
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Figure 27: Contributions to the strength function G T (p) versus the nuclear density p. 



The two- and three-body contributions to the isovector strength functions A(p), G T (p), Gd(p), G so (p) 
and Gj(p) follow the same pattern as outlined in the previous section for the isoscalar strength func- 
tions. The corresponding analytical expressions can be found in ref. [90]. An interesting observation is 
that the Fujita-Mayazawa mechanism for induced spin-orbit couplings from three-body forces is not op- 
erative in the isovector channel, since the 27r-exchange three-body Hartree diagram leads to a vanishing 
contribution to G so (p). 

The results for the isosvector strength functions are shown in Figs. [2614301 The following features 
can be observed. The Hartree-Fock approximation seems to work better for isovector quantities. The 
asymmetry energy at saturation density comes out as A(po) = 26.5 MeV, which agrees within 20% 
with the empirical value (35 ± 2) MeV. In fact, a second-order calculation |106j with two-nucleon forces 
increased the asymmetry energy by approximately 20% and reached the empirical value. The strength 
functions G T (p) and Gy(p), related to the splitting of the proton and neutron effective masses and 
the isosvector surface term (Vp p — Vp n ) 2 , are typically one order of magnitude smaller than their 
isoscalar analogs F T (p) and -Fv(p)- Due to the absence of a large three-body term, the isovector 
spin-orbit coupling strength G so (p) is now dominated by the short-distance contribution (namely, the 
low-energy constant C 5 /8) and thus has only a very weak density dependence. The quadratic spin-orbit 
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Figure 28: Strength function Gy(p) of the isovector surface term (Vp p — Vp n ) 2 versus the nuclear 
density p. 
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Figure 29: Strength function G so (p) of the isovector spin-orbit coupling term (Vp p — Vp r , 
versus the nuclear density p. 
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versus the nuclear density p. 
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coupling strength Gj(p) is again dominated by l7r-exchange contribution which induces a strong density 
dependence. It should be noted that within Skyrme phenomenology the isovector part of the energy 
density functional is presently not well determined. For example, no definite choice could be made 
in ref. |108j between different density dependences (~ p p + 7p n ,7 = 0, 1,2) for the neutron spin-orbit 
potential in Pb isotopes. 

The physics of nuclear structure at large neutron excesses will be explored in the near future at 
rare isotope facilities. The predictions obtained from chiral two- and three-body interactions with their 
definite isospin-dependence can serve as a guideline for the theoretical exploration of this field. 

7 Chiral Effective Interactions and Fermi Liquid Theory 

Landau's theory of normal Fermi liquids [1U9] has long served as a standard basis for understanding 
strongly-interacting Fermi systems at low temperatures. Although more than fifty years old, Landau's 
original insights foreshadowed key ideas in the development of effective field theories |110j and the 
renormalization group |111[I112JI113| . In Fermi liquid theory, the relevant degrees of freedom associated 
with low-energy, long- wavelength excitations of a many-body Fermi system are quasiparticles interacting 
in the vicinity of the Fermi surface through a residual force. Although the quasiparticle interaction is 
not necessary weak, small perturbations of the system excite relatively few quasiparticles, which then 
forms the basis for an expansion of the theory in terms of the quasiparticle density. 

The microscopic foundation for Fermi liquid theory is based on an analysis of the one- and two- 
particle Green's functions |114j . In the vicinity of the Fermi surface, the one-body Green's function 
can be decomposed into quasiparticle and background contributions, and the effects of the background 
contribution are absorbed into the effective coupling between quasiparticles. In the original applications 
of Fermi liquid theory to finite nuclei [115[ 1116] , the quasiparticle couplings were adjusted to empirical 
data, resulting in a phenomenological theory capable of linking seemingly disconnected phenomena. 
The perturbative approach to nuclear Fermi liquid theory based on microscopic models of the strong 
nuclear force was initiated by Brown and collaborators [117] and has in recent years been updated to 
take advantage of the improved convergence properties of chiral and low-momentum nuclear two- and 
three-body forces pH [TEHJ [1061 E2B Q2U • 

In the present section we review the Fermi liquid description of nuclear and neutron matter in the 
context of chiral effective field theory. A microscopic approach framed in many-body perturbation theory 
is shown to provide a successful description of many nuclear matter observables when the leading-order 
medium corrections from two- and three-body forces are included. The extrapolation to pure neutron 
matter is considered, and applications relevant for neutron star structure and evolution are discussed. 
We focus on the convergence of the perturbative expansion with chiral nuclear forces and the role of 
chiral three-nucleon forces. For a more general review of the Fermi liquid approach to nuclear many- 
body systems, we refer the reader to ref. [122J. 

7.1 Symmetric Nuclear Matter 

Numerous properties of isospin-symmetric nuclear matter in the vicinity of the saturation density po 
are well constrained by empirical data and serve as a benchmark for microscopic approaches to nuclear 
structure. Within the framework of Fermi liquid theory, specific nuclear observables are directly related 
to the central components of the quasiparticle interaction, which has the general form 

1 oo 

^cent(pi,p 2 ) = ^Y.[ F L + F' L ri-r 2 + (G L + G' L n • f 2 ) at ■ a 2 ] P L (cos6), (100) 

where N Q = 2M*kf/n 2 is the density of states at the Fermi surface, and we have set |pi| = \p 2 \ = kf, 
which allows one to write the quasiparticle interaction in terms of cos 9 — p\-p 2 . The noncentral parts of 
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Figure 31: Diagrams contributing to the second-order quasiparticle interaction (all interactions repre- 
sented by wavy lines are antisymmetrized): (a) particle-particle diagram, (b) hole-hole diagram, and 
(c) + (d) particle- hole diagrams. 

the quasiparticle interaction, which are particularly relevant for neutron matter, will be discussed in the 
following section. The relationship between individual Fermi liquid parameters and nuclear observables 
are well known |116| : 



Quasiparticle effective mass : 



M* 



M N 
Compression modulus : /C 

Isospin asymmetry energy : 

Orbital g-factor: gi 

Spin-isospin response : g' NN 



1 + Fi/3, 



~M* 

Kp 



1 + ^0 



^6(1 + ^/3) 3 ' 



4M^ 



9, N N 



(101) 



Within the framework of many-body perturbation theory, the quasiparticle interaction is derived 
by functionally differentiating the total energy density twice with respect to the particle occupation 
densities 5n 



p,s,t- 



<5£ = Y1 e P$ n pst + - Yl J r (p 1 s 1 t 1 ;p 2 s 2 t2) Snp 18ltl 6 



pst 



2 .- 

p 1 s 1 t 1 
P2«2*2 



n 



P2S2t2 



+ ■■ 



(102) 



where Sj and ti label the spin and isospin quantum numbers. The first- and second-order contributions 
are shown diagrammatically in Fig. EU As shown in the first row of Table [U the quasiparticle interaction 
computed at leading order with chiral two-body forces shares many features with previous calculations 
employing G-matrix effective interactions. In particular, the quasiparticle effective mass M* jMjq = 0.70 
is significantly less than unity, the nuclear symmetry energy /3 = 21MeV is too small compared to the 
empirical value (3, 



erap 



30 — 35MeV, and the compression modulus /C = — 44MeV is negative and 
implies the instability of nuclear matter to density fluctuations. This last feature was quite common in 
previous calculations of the Fermi liquid parameters with G-matrix effective interactions |123[ I124[ 1125] 
and was remedied with the inclusion of virtual collective excitations, represented diagrammatically as 
the sum of particle-hole bubble diagrams and encoded in the Babu-Brown induced interaction [126J, 
which provided sufficient screening in the scalar-isoscalar channel to render nuclear matter stable. 

Recent chiral effective field theory calculations have systematically studied the second-order con- 
tributions to the quasiparticle interaction, given in rows 3-5 of Table [TJ The particle-particle term 
(represented diagrammatically Fig. I3l7 a)) is found to enhance the attraction in the spin- and isosp in- 
independent interaction, while the particle-hole terms (Fig. I3l7 c) and (d)) are found to be strongly 
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Chiral 2N + 3N interactions at &f = 1.33 fm l 




/o [fm 2 ] 


9o [fm 2 ] 


fo [fm 2 ] 


9o [fm 2 ] 


h [fm 2 ] 


01 [fm 2 ] 


f[ [fm 2 ] 


g[ [fm 2 ] 


yd) 
V 2N 


-1.274 


0.298 


0.200 


0.955 


-1.018 


0.529 


0.230 


0.090 


T/(2-pp) 

V 2N 


-1.461 


0.023 


0.686 


0.255 


0.041 


-0.059 


0.334 


0.254 


x A2-hh) 

V 2N 


-0.271 


0.018 


0.120 


0.041 


0.276 


0.041 


-0.144 


-0.009 


,A2-ph) 
V 2N 


1.642 


-0.057 


0.429 


0.162 


0.889 


-0.143 


0.130 


0.142 


V 3N 


1.218 


0.009 


0.009 


-0.295 


-0.073 


-0.232 


-0.232 


-0.179 



Table 1: Fermi liquid parameters (L = 0, 1) computed at nuclear matter saturation density from the 
Idaho N 3 LO chiral nucleon-nucleon potential as well as from the leading N 2 LO chiral three-nucleon 
force. 



repulsive. Together these two contributions largely cancel in the calculation of the F Landau parame- 
ter, and since the second-order particle-hole diagram accounts for much of the full Babu-Brown induced 
interaction, the stabilization mechanism against density fluctuations for chiral nuclear interactions must 
be qualitatively different than in the case of the G-matrix effective interactions. As seen in Table Q] 
both the /i and /q Landau parameters receive coherent contributions from all second-order diagrams. 
This results in a quasiparticle effective mass at the Fermi surface that is nearly equal to the free-space 
nucleon mass. Taken alone, this would substantially decrease the isospin asymmetry energy. However, 
the second-order contributions also play a dominant role in enhancing the isospin dependence of the 
quasiparticle interaction. At nuclear matter saturation density, the resulting isospin asymmetry energy 
is P = 34MeV, which is in very good agreement with the empirical value. These results obtained from 
the Idaho chiral two-nucleon force are qualitatively the same those from evolved low-momentum NN 
interactions [106]. 

In the chiral effective field theory description of nuclear interactions, three-body forces enter at 
N 2 LO in the chiral power counting. The first-order contributions from the leading-order three-nucleon 
force are shown diagrammatically in Fig. [T2] and have been computed analytically for the first time in 
ref. |120j . The surprising feature revealed in these calculations is the substantial repulsion present in 
the spin- and isospin-independent quasiparticle interaction resulting from two-pion exchange dynamics 
and to a lesser extent the chiral three-nucleon contact interaction. As shown in the left panel of Fig. [321 
the Landau parameter /o grows strongly with the density and provides more than sufficient repulsion 
to stabilize nuclear matter at the empirical saturation density, resulting in a compression modulus of 
/C = MeV. This qualitative behavior is foreseen already in calculations of the nuclear matter equation 
of state with chiral and low-momentum interactions, which require three-nucleon forces in order to 
achieve saturation at the empirical density [99]. However, the magnitude of /o relative to the other 
Landau parameters, as seen in Fig. [32] and compiled at nuclear matter density in the last row of Table 
[U is unexpectedly large. The momentum dependence of the quasiparticle interaction in this channel is 
especially weak, producing only a very small decrease in the quasiparticle effective mass. The scalar- 
isovector and vector-isoscalar interactions are identical and nearly zero when averaged over all angles. 
Three-nucleon forces therefore play a minor role in determining the isospin asymmetry energy of nuclear 
matter, which was already well described with two-nucleon forces alone. 

In the right panel of Fig. [32] we show the density dependence of the full quasiparticle interaction 
for an evolved low-momentum nucleon-nucleon potential combined with a three-nucleon force whose 
low-energy constants are refit at the low-momentum cutoff scale [99] . Qualitatively similar results hold 
for the unevolved Idaho nucleon-nucleon interaction and associated three-nucleon force, but the varia- 
tion in cutoff scale and low-energy constants provides an important means for studying the theoretical 
uncertainties present in the perturbative calculation. From Fig. [32] one sees that the density dependence 
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Figure 32: Left panel: density-dependent L = 0, 1 Fermi liquid parameters obtained from the leading- 
order chiral three- nucleon force with low-energy constants given in ref. [32]. Right panel: density- 
dependent L = 0, 1 Fermi liquid parameters computed from low-momentum chiral two-nucleon and 
three-nucleon interactions at the resolution scale A = 2.1fm _1 . 

of the spin- and isospin-independent quasiparticle interaction is particularly strong. This suggests that 
obtaining a quantitatively correct description of the effective mass and compressibility of nuclear matter 
at the saturation density is nontrivial. In contrast, the remaining components of the quasiparticle inter- 
action have a weak and strikingly similar density dependence. Although most properties of symmetric 
nuclear matter are well reproduced as discussed above, the anomalous orbital ^-factor Sgi ~ 0.07 ±0.02 
is considerably weaker than the value Sgi texp ~ 0.23 ± 0.03 extracted from an experimental analysis of 
the isovector giant dipole sum rule [127j : 



dua(El) 



27raNZ 
M N A 



;i + 25gi 



(103) 



where a = 1/137. Overall, however, the perturbative calculation including the first-order (free-space) 
contribution from two-body forces together with the leading Pauli-blo eking and core polarization effects 
from two- and three-body forces provides a successful description of isospin-symmetric nuclear matter 
in the vicinity of the empirical saturation density. 

7.2 Neutron Matter 

In contrast to the case of isospin-symmetric nuclear matter, the quasiparticle interaction in pure neutron 
matter is largely unconstrained empirically but is highly relevant for neutrino processes in neutron 
stars [1281 11291 11301 1131] as well as the response of neutron star matter to strong magnetic fields 
|132[ 11331 11341 1135] . Previous calculations with microscopic and phenomenological nuclear potentials 
have given qualitatively different predictions for the magnetic susceptibility of neutron matter, with 
some phenomenological potentials even giving rise to a ferromagnetic phase transition at several times 
nuclear matter saturation density |135[ 11361 1137] . Both the magnetic susceptibility of neutron matter 
as well as neutrino absorption, emission, and elastic scattering rates depend on noncentral components 
of the quasiparticle interaction |138j 1139] as discussed in refs. |133U130]I131] . Recently, these questions 
have been revisited within the framework of chiral effective field theory including three-nucleon forces 

cm 

In the long- wavelength limit, the quasiparticle interaction in pure neutron matter has the general 
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-0.500 
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0.159 


-0.364 


0.337 
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h [fm 2 ] 
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0.063 





/ [fm 2 ] 





0.135 


-0.168 





-0.031 


-0.134 





-0.279 


0.083 



Table 2: Fermi liquid parameters for the quasiparticle interaction in neutron matter at a density cor- 
responding to a Fermi momentum of kf = 1.7 fm -1 . The low-energy constants of the N 2 LO chiral 
three-nucleon force are chosen to be c\ = —0.81 GeV -1 and C3 = —3.2 GeV -1 . 



form [139] : 



F{p\,P%) = f{Pi,P2)+g(Pi,P2)vi-V2 + h(p 1 ,p 2 )S 1 2(q) + k(p 1 ,p 2 )S l2 (P) 

+l(pi,P2)(vi x 02) ■ {q x P), 



(104) 



where q = p% — p 2 denotes the momentum transfer in the exchange channel, P = p\ + p 2 is the 
center-of-mass momentum of the quasiparticle pair, and $12(10 defines the tensor operator S\ 2 (v) = 
3<7i • v a 2 ■ v — d\ ■ o 2 . The quasiparticle interaction in eq. fjl04j) is invariant under parity, time- reversal, and 
transformations that interchange particle labels. The presence of the medium allows for the possibility 
of terms that break Galilean invariance and that depend explicitly on the center of mass momentum 
P; namely, Si 2 (P) and Ai 2 (q,P) = (a\ x <j 2 ) ■ (q x P). In ref. |121] a practical method has been 
developed for extracting the various scalar functions from a spin-space decomposition of the quasiparticle 
interaction. Additionally, three-nucleon force contributions to the noncentral quasiparticle interaction 
were computed for the first time. 

The central quasiparticle interaction in neutron matter is qualitatively similar to that described in 
the previous section for isospin-symmetric nuclear matter (with Ti ■ r 2 —¥ 1). The individual contri- 
butions from two- and three-nucleon forces to the L = 0,1,2 Landau parameters at various orders in 
perturbation theory for kf = 1.7fm _1 are shown in Table [2j The first-order contribution from two- 
nucleon forces is simply a kinematically constrained version of the free-space interaction and therefore 
must be independent of the total quasiparticle momentum P. Second-order contributions give rise to all 
possible noncentral interactions, which from Table |2] are seen to be significantly smaller than the central 
components. The exchange tensor interaction present in the free-space nucleon-nucleon interaction is 
largely reduced by Pauli-blocking effects in the nuclear medium. The center-of-mass tensor interaction 
arises solely from Pauli-blocking and core polarization effects at second-order in perturbation theory 
with two-body forces, while the novel spin-nonconserving cross- vector interaction Ai 2 (q, P) is generated 
both from two- and three-body forces, whose momentum-averaged contributions are of comparable 
strength. 

At leading order the compressibility /C of neutron matter at po is unphysically small in comparison 
to auxiliary-field diffusion Monte Carlo calculations with realistic two- and three-nucleon forces which 
find that at p = 0.16 fm -3 , /C ~ 520 MeV [136]. Second-order contributions from chiral two-body 
forces do little to change this picture, but as in the case of symmetric nuclear matter the leading-order 
chiral three-nucleon force provides substantial additional repulsion, yielding a compression modulus 
/C = 550 MeV. The f\ Landau parameter governing the quasiparticle effective mass at the Fermi surface 
again receives a nearly negligible contribution from three-body forces, but strong polarization effects 
at second order increases the effective mass from M*/Mjy = 0.82 to M*/Mjy = 1.04. As seen in Fig. 
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Figure 33: Energy per particle of neutron matter from both chiral and low-momentum two- and three- 
body interactions. The cutoff scale associated with the low-momentum interaction is A = 2.1 fm -1 . 
The curve labeled 'APR' is taken from the variational calculations of Akmal et al. [54\ , and the curve 
labeled 'QMC is taken from the quantum Monte Carlo calculations in refs. [1401 14"9] . 



1341 the Landau parameters fi associated with the spin-independent quasiparticle interaction depend 
sensitively on the neutron density. 

Neglecting effects from noncentral interactions, the magnetic susceptibility is determined by the 
Landau parameter G : 

o N 



X = Vn 



1 + Gn 



(105) 



where //„ = —1.913 is the free-space neutron magnetic moment in units of the nuclear magneton, and 
Go — ^ofl'o- Noncentral components in the quasiparticle interaction result in effective charges (magnetic 
moments) that are not scalars under rotations of the quasiparticle momentum. The spin susceptibility 
then involves separately the longitudinal and transverse components of the magnetic moment [132, 133]. 
The Landau parameter go decreases with the nuclear density as shown in Fig.[3U but there is no evidence 
for a spin-instability in the vicinity of nuclear matter saturation density. 

As a consistency check on the perturbative calculation of the quasiparticle interaction in neutron 
matter, we reproduce in Fig. [33] the neutron matter equation of state at zero temperature computed to 
the same order in many-body perturbation theory. The results agree well with variational and Monte 
Carlo calculations of neutron matter employing realistic and phenomenological nuclear forces [Ml I140J, 
H9"] . The recent neutron matter calculation [69J including the N 3 LO chiral three- and four-neutron forces 
gives results that are very similar to those shown in Fig. |33j The dependence on the resolution scale 
arises mostly from the different choice of the low-energy constant C3 in the two calculations. Recently, 
it has been shown that refitting C3 to peripheral partial waves with the running of the cutoff scale yields 
a much reduced regulator dependence of the neutron matter equation of state [70J . 

The dynamics of supernovae explosions are particularly sensitive to neutrino processes occurring in 
the neutrinosphere, a warm and low-density (n ~ po/10) gas of neutron-rich matter. At such small 
densities, three-nucleon correlations are expected to be negligible, and an accurate understanding of 
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Figure 34: (Color online) Density-dependent Fermi liquid parameters inlcuding first- and second-order 
contributions from the chiral N 3 LO nucleon-nucleon potential of ref. [31] as well as the N 2 LO chiral 
three-nucleon force to leading order. 
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Figure 35: The rms energy transfer in neutrino scattering from neutron matter as a function of the 
temperature, density, and neutrino energy. Results are computed with chiral two-body forces at different 
orders in the chiral power counting. Figure reproduced from ref. [131] . 
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neutrino scattering, production, and annihilation rates is accessible from two-nucleon forces alone. 
These processes determine the supernova neutrino spectrum and may be important for heating effects 
that re-energize the shock wave. The typical scale for neutrino energies e v ~ 30 MeV and momenta p v ~ 
0.15 fm _1 justifies a treatment within linear response theory. Supernova simulations generally include 
neutrino-pair bremsstrahlung and absorption rates computed in the one-pion exchange approximation. 
For degenerate conditions, where the temperature is small relative to the neutron Fermi energy, Landau's 
Fermi liquid theory is applicable has been used in recent work |131j to compute spin-spin response 
function of neutron matter employing chiral two-body interactions at N 3 LO. 

Within Fermi liquid theory, the density-density and spin-spin response functions are computed by 
solving the transport equation for quasiparticles. The relaxation rate that is used to approximate the 
collision integral in the Landau transport equation involves contributions from central and noncen- 
tral components of the quasiparticle interaction, the latter being essential for describing two-nucleon 
response. Contributions from chiral nuclear interactions beyond the leading-order one-pion exchange 
potential substantially reduce the relaxation rate, which consequently reduces the mean-square energy 
loss in neutrino scattering, as shown in Fig. 



8 Nuclear Chiral Thermodynamics 

8.1 Nuclear Phase Diagram and Liquid-Gas Transition 

The thermodynamic properties of nuclear and neutron matter play an important role in applications 
to heavy-ion collisions and astrophysics. The plateau observed in the caloric curve of the nuclear 
fragments in nucleus-nucleus collisions can be viewed as the trace of a first-order liquid-gas phase 
transition |141[ 11421 1143] . In the context of nuclear astrophysics, the recent observation of a two-solar- 
mass neutron star |144] strongly constrains the equation of state of hadronic matter and rules out 
various exotic models which tend to produce softer equations of state. 

Previous chapters have introduced and developed the framework of chiral effective field theory 
and its application to the nuclear many-body problem at zero temperature. In the present section 
we review progress that has been achieved in extending this framework to finite temperatures. The 
previous ordering scheme for the nuclear energy density in terms of the number of medium insertions 
is now generalized and applied to the free energy density as a function of density p and temperature T. 
Introducing the free energy per nucleon, F(p,T), the free energy density can be written as the sum of 
convolution integrals: 

oo oo oo 

pF(p,T) = 4jdppK,i(p)n(p) + jdpijdp2)C2(pi,P2)n(p 1 )n(p2) 

o oo 

oo oo oo 

+ J dpi J dp 2 j dp 3 K. 3 (p 1 ,p2,P3)n(p 1 )n(p 2 )n(p 3 ) + pA(p,T) , (106) 



where 



n(p) 



p ., ^p 2 /2M N - p]- 1 



2vr 2 



1 + exp ■ 



T 



(107) 



is the density of nucleon states in momentum space. It is the product of the temperature- dependent 
Fermi-Dirac distribution and a kinematical prefactor p/2n 2 that has been included in n(p) for conve- 
nience. Complete expressions of the kernels /C, can be found in |145j . The term A(p,T) is called the 
anomalous contribution. It is associated with a Fock term involving second order l7r-exchange [146] and 
arises from the smoothing of the Fermi surface at nonzero temperature. In the range of temperatures 
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considered here the effect of the anomalous term is small. The one-body effective "chemical potential" 
\x in the distribution n{p) is related to the density by 



p = 4 / dppn(p) . (108) 



The kernel /Ci is the non-interacting nucleon gas contribution, while the kernels /C2 and /C3 encode 
the effects of the interactions among two and three nucleons. The pressure is then computed from the 
standard thermodynamic relation 



2 dF(p,T) 



P(p,T) = p* f p '- . (109) 

In addition to isospin-symmetric nuclear matter, the extension to asymmetric nuclear matter and 
neutron matter has also been explored. The only changes required involve isospin factors and the 
introduction of separate proton and neutron thermal occupation probabilities d p (pj) and d n (pj). Note 
that for pure neutron matter the Pauli principle forbids a (momentum independent) three-body contact 
interaction. The behavior at very low density is governed by the large neutron-neutron scattering 
length, a nn ~ 19 fm [IE], and as discussed in Section [5.21 an all-order resummation is required for a 
realistic description at low densities (p„ < 0.02 fm -3 ). 

In this section we first concentrate on the equation of state of isospin-symmetric nuclear matter for 
different temperatures, featuring the liquid-gas phase transition. The calculations are then extended to 
isospin-asymmetric matter, focusing on the change of thermodynamic properties with varying proton- 
neutron asymmetry. The topics of interest are the asymmetry free energy, its dependence on density 
and temperature and the validity of the (empirically observed) parabolic approximation for the free 
energy as a function of the asymmetry parameter, 5 = (p n — p p ) /(p n + p p ). 

The free energy per particle F(p, T) is calculated using eq. (11061) as a function of density p for 
a sequence of temperatures up to 25 MeV, with the input for the interaction kernels K n specified in 
ref. |145| . The result is shown in Fig. |36j The dotted lines indicate the non-physical behavior of the 
equation of state in the first-order liquid-gas transition region. This unphysical part is then replaced 
by the physical one (solid lines) obtained from the Maxwell construction. At zero temperature the free 
energy is equal to the internal energy of the system, which has the well-known saturation point located 
at Eq = —16.0 MeV and po = 0.16 fm -3 , corresponding to a Fermi momentum of fc/o = 1.33 fm - = 262 
MeV. At finite temperatures the free energy develops a singular behavior for p — y 0, a feature that is 
present in numerous other many-body calculations [HH 11471 . 



In Fig. [33 we show isotherms of the pressure P(p, T) as a function of the temperature and density. 
The picture that emerges is reminiscent of a Van der Waals gas with its generic first-order liquid-gas 
phase transition. As outlined in Section 4.3, chiral nucleon-nucleon dynamics generates intermediate- 
range attractive interactions (e.g., from 27r-exchange with intermediate A-excitations) that exhibit a 
characteristic e~ 2m,Tr /r 6 behavior at distances on the order of 1 — 2fm (see ref. |40j). Such mechanisms 
provide nearly half of the attraction required to bind nuclear matter at zero temperature. The liquid-gas 
phase transition then results from a sensitive balance between intermediate range attraction and short- 
range repulsion, the latter encoded in contact terms representing dynamics unresolved at the scales 
relevant for typical nuclear Fermi momenta kf considered here. 

The critical temperature for the liquid-gas phase transition is observed at T c ~ 15.1 MeV. For 
temperatures less than the critical temperature T c , the usual Maxwell construction is applied, where 
the pressure is kept constant in the liquid-gas coexistence region. Empirical estimates of the critical 
temperature are deduced from fission and multi-fragmentation measurements, which indicate a value 
of T c ~ 15 — 20 MeV |148j . Similar values for the critical temperature are obtained from calculations 
|149] employing phenomenological Skyrme interactions. 
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Figure 36: Free energy per particle of isospin-symmetric nuclear matter as a function of the temperature 
and baryon density p. The dotted line denotes the non-physical part of the free energy in the liquid-gas 
coexistence region. The physical free energy (solid line) is obtained using the Maxwell construction. 
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Figure 37: Pressure isotherms as a function of density and temperature for isospin-symmetric nuclear 
matter displaying a first-order liquid-gas phase transition. The dotted lines at low temperature indi- 
cate the non-physical behavior of the isotherms in the transition region, and the physical pressure is 
calculated from the Maxwell construction. The dashed line delineates the boundary of the coexistence 
region, and the dot denotes the critical point at T c ~ 15.1 MeV and p c ~ po/3. 
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In Fig. |38]we show the dependence of the pressure on the nucleon chemical potential (including the 
free nucleon mass), 



/! = M N + 



^p|)f(p,d 



(110) 



The non-physical part of the pressure is shown as the dotted curves in Fig. [37] and manifested in 
the double-valued behavior of P at temperatures below the critical temperature. In the coexistence 
region, the physical pressure and chemical potential are constant and given according to the Maxwell 
construction. The temperature at which the pressure becomes single-valued is the critical temperature 
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Figure 38: Dependence of the pressure isotherms on the nucleon chemical potential p for symmetric 
nuclear matter. The non-physical behavior of the equation of state in the liquid-gas coexistence region 
for temperatures below T c is denoted with the dotted lines. The physical pressure and chemical potential 
in the coexistence region are constant and determined from the Maxwell construction. The critical point 
is denoted with the large dot. 




The T — p phase diagram, shown in Fig. [39], collects the relevant information regarding the liquid-gas 
phase transition. The first-order transition region ends at the critical point (T c = 15.1 MeV) and is 
denoted with the dot. The associated critical values of the pressure, baryon chemical potential, and 
density are found to be P c ~ 0.261 MeVfm -3 , p c ~ 914.7 MeV, and p c ~ 0.053 fm" 3 . The T- p phase 
diagram indicates that the liquid-gas coexistence region extends over a wide range of densities up until 
the Fermi liquid phase is realized at po — 0.16 fm~ 3 , the equilibrium density of nuclear matter. 

At zero temperature the third law of thermodynamics provides important constraints on features of 
the T — p diagram. In particular, at T = the boundary of the phase coexistence region has an infinite 
slope in the T — p diagram. Additionally, the chemical potential at zero temperature is given by the 
total energy per particle at the saturation point, that is, p = M^ + E ~ 923 MeV. 



8.2 Is o spin- A symmetric Nuclear Matter 

It is now of interest to investigate how the liquid and gas phases characteristic of isospin-symmetric 
(N = Z) nuclear matter evolve with changing proton fraction, x p = Z/A. The properties of isospin- 
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Figure 39: The phase diagram of isospin-symmetric nuclear matter. The critical point is indicated by 
the dot on the boundary of the phase coexistence region. 



asymmetric matter with 0.5 > x p > are governed by the detailed isospin dependence of the chiral NN 
interaction. It is almost entirely controlled by the isospin dependence of one- and two-pion exchange 
processes once the relevant nucleon-nucleon contact terms are fit to the empirical isospin asymmetry 
energy A(p ,T = 0) ~ 34MeV. 

Consider the evolution of the saturation point, defined as the minimum of the energy per nucleon 
at T = 0, as the proton fraction decreases. The result is shown in Fig. HUJ Starting from its value 
for symmetric nuclear matter, E ~ —16 MeV at p — 0.157 fm -3 , the binding energy per particle 
continuously reduces in magnitude with decreasing x p until it vanishes for a value of the proton fraction 
x p ~ 0.12, beyond which the neutron-rich matter is unbound at zero temperature. 

In Fig. HTJwe show the phase diagram as a function of the proton fraction x p . The dashed line shows 
the trajectory of the critical point and its disappearance for proton fractions below x p ~ 0.05. At this 
value of the proton fraction, the coexistence region reduces to a single point, indicating the absence of 
a liquid-gas phase transition. Neutron-rich matter with x p < 0.05 therefore exists only in a gaseous 
phase. We note as well that the critical point disappears at a small but nonzero value of the pressure. 

From Fig. HT]we see that for proton fractions above x p ~ 0.1, the liquid-gas coexistence region starts 
at p = 0. At T = and from x p — 0.1 to 0.5, the boundary between the coexistence and Fermi liquid 
regimes covers a density range between about 0.5 p and nuclear matter saturation density density p . 
The behavior of the coexistence region changes qualitatively for a proton fraction x p = 0.12, which is 



the point at which the zero-temperature binding energy vanishes. From Fig. HO] we observe that for 
x p < 0.12 the absolute minimum of the energy per particle is located at p = but there may still exist 
a local minimum in E(p,x p ) at a finite value of the density. This indicates that neutron-rich nuclear 
matter with values of x p < 0.12 is in the gaseous phase at very low density and T = 0. As the density 
increases, the neutron-rich matter then enters the coexistence region. For values of the proton fraction 
0.05 < x p < 0.12, nuclear matter will not be self-bound but may nevertheless still have a liquid-gas 
phase transition. 

The treatment of nuclear matter chiral thermodynamics presented thus far is oversimplified at low 
densities where nuclear clustering occurs (mainly driven by the Coulomb repulsion of the protons). Pre- 
vious work |150] that combines the appearance of deuteron, triton and helium clusters with relativistic 
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Figure 40: Dependence of the energy per nucleon (and saturation point) of nuclear matter on the proton 
fraction x p at T = 0. The dashed line denotes the evolution of the saturation point with varying proton 
fraction. For proton fractions below x p = 0.12, the energy remains positive at all densities. 
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Figure 41: Dependence of the nuclear matter phase diagram on the proton fraction x p . The evolution 
of the critical point is given by the dashed line. 
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mean-field theory suggests, however, that only modest changes in the T — p phase diagram arise. When 
clustering effects are included, the position of the critical point changes by less than 1% in p c and less 
than 10% in T c . 

The free energy per particle can be expanded in powers of the nuclear asymmetry parameter 5 = 
(pn — Pp)/p — 1 ~ ' 2a; p around the free energy at 5 = 0: 

F(p p , p n , T) = F(p, T) + A(p, T) 6 2 + 0(5 A ) , (111) 

which defines the asymmetry free energy per particle A(p,T). If isospin-symmetry breaking effects 
are ignored, the expansion of the free energy per particle includes only even powers of 5. In this 
approximation, nuclear matter is invariant under the interchange of neutrons and protons. 

The validity of the parabolic approximation has been shown in ref. |145j . where the free energy 
difference relative to isospin-symmetric nuclear matter as a function of 5 2 for different nucleon densities 
P = Pn + Pp has been analyzed. At T = the quadratic dependence on 5 holds very well even up to large 
values of the isospin asymmetry parameter. For higher temperature (e.g., T = 20MeV) one observes a 
slight bending, especially at low density. 

At nuclear matter saturation density po, the contact terms of the isospin-dependent part of the 
nucleon-nucleon interaction have been fit to the isospin asymmetry energy A(p Q ,T = 0) ~ 34.0 Me V. 
Previous determinations from fits of nuclear masses (131 1151] as well as relativistic mean-field models 
constrained by the properties of specific nuclei |152] or by giant dipole resonances |153] give values of 
A(p ) in the range 33 — 37MeV. The asymmetry energy at a lower density, p = 0.1 fm -3 , has also been 
estimated in the latter paper and yields a value for the isospin asymmetry energy between 21 and 23 
MeV. This range is slightly below the value A(p = 0.1fm _3 ,T = 0) ~ 23.9 MeV computed within the 
present chiral effective field theory framework. 

In the vicinity of the nuclear matter saturation density po, the asymmetry energy at zero temperature 
can be expanded in powers of p — po : 

„ „ ts /„ „ \ 2 

(112) 

The coefficients of the linear and quadratic terms (in p — p$) are found in the present case to be L ~ 90 
MeV and K as ~ 153 MeV. In particular, this value of L is consistent with empirical constraints from 
isospin diffusion, which give L = 88 ± 25 MeV [154j . 

For small values of the isospin asymmetry 5, the saturation density is reduced to po(l — 3L5 2 /K). 
The nuclear compression modulus K(8) is then often expressed in terms of an expansion in powers of 
5: 

K(8) =K + K T 5 2 + C(5 4 ) , K T = K as -QL , (113) 

where K is the compression modulus of isosopin-symmetric nuclear matter, and K T is referred to 
as the isobaric compressibility. The value calculated in the present framework is K T = —388 MeV. 
Recent measurements of giant monopole resonances in even-mass-number isotopes yields the result 
K T = -550 ± 100 MeV [T55] . 

8.3 Comparisons to other Approaches 

In the present section we discuss recent calculations of the thermodynamic properties of nuclear mat- 
ter within the self-consistent Green's function method or mean-field theory |1561 11571 1158] 1159] [160J. 
As noted in ref. |156] , many previous studies of nuclear matter at nonzero temperature carried out in 
the mean-field approximation treat the temperature dependence in a simplistic way, namely, the zero- 
temperature step-function momentum distribution is just replaced by the corresponding Fermi-Dirac 
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Figure 42: Energy per particle of neutron matter calculated at temperatures T = 3, 6 and 10 MeV as 
function of density, using low momentum NN and 3N interactions. Left: Hartree-Fock results. Right: 
results including 2nd order NN interactions. Also shown for comparison are curves from a many-body 
calculation (FP) |162] and a model-independent virial expansion (virial) |147] at low density. 

distribution. Such an approach neglects the temperature dependence of the phenomenological inter- 
actions that account for nuclear correlations. In microscopic many-body methods, on the other hand, 
Pauli-blocking effects generated by the nuclear medium are weakened as the temperature increases. 
Hence, nuclear matter properties and correlations depend on the temperature in a non-trivial way. 

In the SCGF formalism, the set of equations relating the in-medium T-matrix, the nucleon self- 
energy, and the single-particle Green's function are solved self-consistently. The equation of state of 
isospin-symmetric nuclear matter has been calculated in this approach |158] starting from the high- 
precision CD-Bonn and Nijmegen nucleon-nucleon potentials. The contributions from three-nucleon 
forces are included through an effective (medium-dependent) two-body interaction. Three-body corre- 
lations are found to strongly affect the pressure P(p,T,S) and result in a liquid-gas coexistence region 
that is much reduced in size. The critical temperatures associated with the CD-Bonn nucleon-nucleon 
potential and the Nijmegen potential are T c = 12.5 MeV and T c = 11.5 MeV, respectively. The critical 
density is found to be p c ~ 0.09 — 0.11 fm~ 3 , while the critical pressure has the value 0.15 MeVfm -3 . The 
results obtained within a chiral effective field theory treatment and described in previous sections are 
therefore significantly different than the results based on the self-consistent Green's function approach 
with one-boson exchange nuclear interactions. 

The properties of nuclear matter have also been calculated |160] within the framework of relativistic 
mean field theory incorporating density-dependent meson-nucleon couplings that account for effects of 
the medium. Such a method has been successful in previous calculations of finite nuclei and infinite 
nuclear matter. The critical temperature T c = 13.2 MeV is found for symmetric nuclear matter, and for 
proton fractions below x p ~ 0.07 the liquid-gas coexistence region vanishes. The picture that emerges 
in these calculations is closer to the one reported in the present review. The main difference between 
the two is that mean field theory gives a slightly smaller phase transition region. 

To close this section we mention the work of ref. |161] in which the energy per particle of neutron 
matter has been calculated at several temperatures T < 10 MeV as a function of density, using a low- 
momentum neutron-neutron interaction that is close in spirit to the chiral interaction discussed in this 
review. Two- and three-body interactions are taken into account. Representative results from these 
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calculations are displayed in Fig. H2J The left panel shows first order (Hartree-Fock) results, including 
a comparison with a model-independent virial expansion at low density. The right panel demonstrates 
significant improvements when the NN interaction is treated in second order. 

9 Thermodynamics of the Chiral Condensate 

The chiral condensate (qq) is defined as the expectation value of the scalar quark density and plays a 
fundamental role in QCD as an order parameter for the spontaneous breaking of chiral symmetry that 
occurs at low energy scales. The dependence of (qq) on the temperature and baryon density is key 
to locating the chiral restoration transition boundary in the phase diagram of QCD. The condensate 
is expected to vanish at high temperatures and/or densities, signaling the crossover from the Nambu- 
Goldstone phase, in which chiral symmetry is spontaneously broken, to the Wigner-Weyl realization of 
chiral symmetry in QCD. 

It is therefore of principal interest to perform a systematic calculation of the thermodynamics of 
the chiral condensate [163J. For such a calculation it is necessary to understand the dependence of the 
free energy density on the light quark mass (or equivalently, the pion mass). In-medium chiral effective 
field theory provides the appropriate framework for such a study, given that the theory provides explicit 
access to the pion mass dependence through one- and two-pion exchange dynamics and the resulting 
two- and three-body correlations in the nuclear medium. 

The starting point is the free energy density, J-"(p, T) = pF(p, T), of spin-saturated isospin-symmetric 
nuclear matter, given by the series of convolution integrals in eq. (I106p . This series is written in terms 
of one-body, two-body and three-body kernels, /Q, /C2 and /C3, respectively. 

The one-body kernel /Ci in eq. f)106p represents the contribution of the non-interacting nucleon gas 
to the free energy density and it reads |146| : 

Mri = M N + „-^-^. (114) 

The first term in /Ci gives the leading contribution (density times nucleon rest mass) to the free energy 
density. The remaining terms account for (relativistically improved) kinetic energy corrections. 

The two- and three-body kernels, /C2 and /C3, specifying all one- and two-pion exchange processes up 
to three loop order for the free energy density, are given in their explicit form in refs. |146[ 139] , 1145] . Let 
us recall that after fixing a few contact terms the free energy density computed from these interaction 
kernels provides a realistic nuclear equation of state up to densities p < 2 p . What matters in the present 
context is the dependence of the kernels /C2 and /C3 on the light quark mass, m q , or equivalently, on the 
pion mass, m n , that is introduced by pion propagators and by pion loops. 

The Feynman-Hellmann theorem is used to establish an exact relation between the temperature- 
and density-dependent quark condensate (qq) (p, T) and the derivative of the free energy density of 
nuclear matter with respect to the light quark mass m q . From the Gell-Mann-Oakes-Renner relation 
m^f% = —m q (0\qq\0), the ratio of the in-medium to vacuum quark condensate is given by 

(qq)(p,T) P dF(p,T) 

(0\qq\0) ~ ft dml ' l ] 

where the derivative of the free energy density with respect to m^ is taken at fixed T and p. Both 
the vacuum quark condensate (0|gg|0) and the pion decay constant f^ are taken in the chiral limit. 
Similarly, m^. represents the leading linear term in the expansion of the squared pion mass in terms of 
the quark mass. 
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The quark mass dependence in the one-body kernel K\ is implicit through its dependence on the 
nucleon mass M N . The condition dp/dM N = leads to the dependence of the effective one-body 
chemical potential jl on the nucleon mass M^: 

^ - 3p fig = -AM N r dp^ . (116) 



dM N 2M N tt'^ u Jo ^ p 

The nucleon sigma term a^ = (N\m q qq\N) = m^dM^ /dml relates the variation of the nucleon mass 
to that of the pion mass. Combining the above two relationships then leads to the expression for the 
derivative of the one-body kernel with respect to the pion mass: 



gjCi _ a N f 3p P 2 gg 

dm\ " m\ I 2M N Q ' 3M N 8M N 



+ — h > (in) 



In the zero temperature limit, the terms in eq. (111T)) provide the linear decrease of the chiral condensate 
with the nucleon density. The corrections from the kinetic energy account for the small difference 
between the scalar and vector densities. The nucleon sigma term is determined empirically at the 
physical pion mass to be &n = (45 ± 8) MeV [26J. Recently, smaller values of the nucleon sigma term 
have been suggested based on lattice QCD calculations of the quark mass dependence of baryon masses 
and accurate chiral extrapolations |165j . The value of a^ determined in this way is consistent with the 
empirical value to within the experimental uncertainty. However, a considerably larger nucleon sigma 
term, <jn — 60 MeV, has recently been discussed in ref. [166] . In the present discussion based on ref. 
|163j . the central value aN = 45 MeV has been employed. 

There are some three-loop contributions that are of special relevance for the in-medium quark 
condensate. The corresponding two-pion exchange mechanisms arise from the chiral symmetry breaking 
7T7tNN contact- vertex proportional to Cim 2 . Note that this parameter is equivalent to the leading 
contribution (linear in the quark mass) to the nucleon sigma term aN = — 4cim 2 + 0(m^). Concerning 
the free energy density pF(p,T) or the equation of state of nuclear matter the c\ml contributions 
are almost negligible. However, when taking the derivative with respect to m^ as required for the 
calculation of the in-medium condensate, these contributions turn out to be of similar importance as 
other interaction terms. The pertinent two- and three-body kernels /C23 can be found in ref. |163j . 

It is meaningful to incorporate also effects from thermal pions. The pressure (or free energy density) 
of thermal pions gives rise, through its m^-derivative, to a further reduction of the T-dependent in- 
medium condensate. In the two-loop approximation of chiral perturbation theory one finds the following 
shift of the chiral condensate ratio due to the presence of a pionic heat bath |167[ I168[ 1169] : 

5(qq)(T) 3ml /m ir \ f m 



^){^4^N^)-^© + ^]} 



(0\qq\0) (^U) 2 K Tn Stt^ 2 

(118) 

where the functions Hi t3 (m n /T) are defined by integrals over the Bose distribution of thermal pions: 



Hi(v)= dx 2 - H 3 (y)=y- 2 dx V - V (119) 

Jy y/x z — y z {e x — 1) Jy e x — 1 

It is worth emphasizing that in-medium chiral perturbation theory with this dynamical input yields a 
realistic nuclear matter equation of state [39 |ll45j . including a proper description of the liquid-gas phase 
transition with critical temperature T c ~ 15 MeV. The three-loop calculation of the free energy density 
is expected to be reliable up to about twice nuclear matter saturation density and for temperatures 
below about T ~ 100 MeV, where the hot and dense hadronic matter remains in the spontaneously 
broken chiral symmetry regime. 
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Figure 43: Density dependence of the chiral condensate at temperature T = 100 MeV. Starting from 
the linear density dependence (dashed curve) associated with the free nucleon Fermi gas, the following 
interaction contributions are successively added: one-pion exchange Fock contribution, iterated one-pion 
exchange, irreducible two-pion exchange, two-pion exchange with intermediate A excitations, two-pion 
exchange contribution proportional to the low-energy constant c\. For the sake of clarity in presenting 
the pattern of medium modifications, the contribution from thermal pions is omitted here 



In Fig. H3] we show a typical example (T = 100 MeV) of the effects of interaction contributions 
to the density dependence of (qq) (p, T) arising from the pion mass derivative of the chiral two- and 
three-body kernels /C2 and /C3. Of particular importance is the pion- mass dependence of the terms 
involving virtual A(1232) excitations, which delay the tendency towards chiral symmetry restoration 
as the density increases. Summing together all one- and two-pion exchange processes contributing to 
d^/drri^ and dlC^/dm^, the chiral condensate at T = 100 MeV is not far from the linear density 
dependence characteristic of a free Fermi gas. However, this behavior is the result of a subtle balance 
between attractive and repulsive correlations and their detailed pion-mass dependences. Iterated one- 
pion exchange and irreducible two-pion exchange alone would have resulted in the system becoming 
unstable not far above normal nuclear matter saturation density as can be seen in Fig. H3j In the chiral 
limit (m,r — > 0) this instability would have appeared even at much lower densities. This emphasizes 
once again the importance of including A(1232)-excitations, and it also highlights the significance of 
explicit chiral symmetry breaking in governing nuclear scales. 

In Fig. HH we show the variation of the chiral condensate with both the baryon density p and the 
temperature T. All nuclear correlation effects are included as well as the small additional shift from 
thermal pions, which is visible only at the highest temperature considered here (T = 100 MeV). The 
actual crossover transition can be extracted from lattice QCD simulations and has the value T ~ 170 
MeV [1701 . 

At zero temperature, the decrease in the chiral condensate beyond nuclear matter saturation density 
is hindered as a result of three-body correlations encoded in /C3 which grow faster than /C2 as the density 
increases. As the temperature rises, the influence of /C3 relative to /C2 is reduced, so that at T = 100 
MeV only a small net effect remains in comparison to the free Fermi gas. 

At T = 0, the solid line in Fig. HH does not yet account for the coexistence region of the nuclear 
liquid and gas phases |145j . Any first-order phase transition is expected to be visible also in other order 
parameters, and the dashed line in Fig. HH based on the usual Maxwell construction indicates this effect. 
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Figure 44: Ratio of the in-medium chiral condensate to its vacuum value as a function of the baryon 
density p and temperature in isospin-symmetric nuclear matter. Effects of thermal pions are included, 
and the dashed line at T ~ results from the Maxwell construction in the nuclear liquid-gas coexistence 
region. 

This feature becomes much more pronounced when the chiral condensate is plotted as a function of the 
baryon chemical potential at low temperatures: 



p = M N +(l + p^-)F(p,T). 



(120) 



The discontinuity associated with the first-order liquid-gas transition at T smaller T c ~ 15 MeV is 
clearly visible in Fig. H5j Another effect generated by the first-order liquid-gas phase transition is that 
the frequently advocated "low-density theorem" requires modification: 



(gg)(p) 

(0|gg|0) 



m lfn' 



121) 



In this equation the effective nucleon sigma term a^ — 36 MeV represents the quark mass dependence 
of the sum M^ + E , where E Q ~ — 16 MeV is the binding energy per particle of nuclear matter at the 
saturation density. The standard version of eq. ( 1 1 2 1 [) contains the nucleon sigma term <jn in vacuum 
and assumes that at sufficiently low densities nuclear matter can be treated as a non-interacting Fermi 
gas. 

The results discussed here set important nuclear physics constraints for the QCD equation of state at 
baryon densities and temperatures that are of interest in relativistic heavy-ion collisions. In particular, 
there is no indication for chiral symmetry restoration at temperatures T < 100 MeV and baryon densities 
at least up to about twice nuclear matter saturation density. The effects of intermediate A(1232)-isobar 
excitations (i.e. the strong spin-isospin polarizability of the nucleon) together with Pauli blocking in 
the nuclear medium play a crucial role in stabilizing the scalar quark condensate at and beyond nuclear 
matter saturation density. 
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Figure 45: Ratio of chiral condensate relative to its vacuum value as function of baryon chemical po- 
tential in isospin-symmetric nuclear matter at low temperatures within the liquid-gas phase coexistence 
region. 

10 Concluding Remarks and Outlook 

Chiral Effective Field Theory, as a low-energy realization of QCD in its sector with light (u and d) quarks, 
is not only a successful framework for the description of pion-pion and pion-nucleon interactions. With 
its systematic counting scheme in powers of "small" Q/A x (where Q stands generically for momentum, 
energy or the pion mass, and A x ~ 1 GeV is the spontaneous chiral symmetry breaking scale), chiral 
EFT has also become the basis for the modern theory of nuclear forces. The present review goes still 
one step further and reports on recent developments concerning the application of chiral EFT to nuclear 
many-body problems, both at zero and finite temperatures. 

The additional "small" scale introduced by the many-body system is the Fermi momentum. Its 
typical magnitude in nuclear physics, kf ~ 2m n , implies that the chiral dynamics of pions in a low- 
energy, low-temperature nuclear environment should also be controlled by an expansion in powers 
of kf/A x . Indeed, the typical average distance between two nucleons in equilibrium nuclear matter, 
d c± 1.8 fm, is again comparable to a "pionic" scale, just slightly larger than the Compton wavelength 
of the pion. Furthermore, another characteristic parameter of nuclear matter, its compression modulus, 
happens to be of a magnitude comparable to twice the pion mass. The only quantity that falls out of this 
pionic scales scenario is the binding energy per nucleon. Its value B/A ~ 16 MeV is abnormally small 
and indicates a quantitative fine-tuning: the well-known subtle balance between intermediate-range 
attraction and strong short-range repulsion in the nucleon-nucleon interaction. 

Details of the short-range NN interaction remain unresolved at momentum scales p < kf. Chiral 
EFT assigns a limited number of NN contact terms to the physics at short distance (r < 1 fm), with 
prefactors (a set of low-energy constants) to be determined by experiment. While the strengths of 
these contact terms generally require non-perturbative resummation techniques to be performed at low 
densities, the intermediate and long range forces governed by one- and two-pion exchange processes in 
the medium are supposed to be accessible to perturbative methods. This is the essence of in-medium 
chiral perturbation theory. In this perturbative hierarchy, three-nucleon forces enter naturally at a 
well-defined order of the low-momentum expansion. 
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Such a conceptual framework turns out to be remarkably successful. Much of the well established 
phenomenology associated with the nuclear many-body problem can now be given more systematic 
foundation that relies substantially on the well-established spontaneous chiral symmetry breaking pat- 
tern of low-energy QCD. The nuclear mean field experienced by nucleons as quasiparticles, including 
its extension to the optical potential for scattering, the density functional approach to finite systems, 
Landau-Migdal Fermi liquid theory, thermodynamics and the phase diagram of nuclear matter featur- 
ing the liquid-gas phase transition - all these important aspects of nuclear theory are covered in this 
report. As a bonus of this framework that permits investigating the energy density and related quanti- 
ties systematically as a function of the pion mass, reliable statements about the behavior of the quark 
condensate at temperatures up to about 100 MeV and densities reaching twice the density of normal 
nuclear matter can be made, indicating stability of this chiral order parameter considerably beyond 
previous expectations that were based merely on the Fermi gas limit. 

Of course, many further reaching questions still need to be addressed. So far, the calculations 
of the energy density (and of the free energy density at finite temperature) have been carried out 
up to three-loop order. This includes important three-body forces but does not (yet) cover issues 
of convergence beyond this order, e.g. concerning the role of four-nucleon correlations. At this point 
Monte Carlo studies of the kind performed in ref. |172] may prove useful in the near future. As a further 
important topic, the quest for a sufficiently stiff equation-of-state satisfying the constraints imposed by 
the existence of a two-solar-mass neutron star |144] is under investigation. 
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